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GIANT MOLECULAR CLOUDS AND STAR FORMATION IN THE NON-GRAND DESIGN SPIRAL GALAXY NGC 6946 
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ABSTRACT 

We present high spatial resolution observations of Giant Molecular Clouds (GMCs) in the eastern part of the 
nearby spiral galaxy NGC 6946 obtained with the Combined Array for Research in Millimeter-wave Astronomy 
(CARMA). We have observed CO(l -> 0), CO(2 ^ 1) and 13 CO(l -> 0), achieving spatial resolutions of 5. "4 x 
5."0, 2 "5 x 2."0 and 5 "6 x 5. "4 respectively over a region of 6 x 6 kpc. This region extends from 1.5 kpc to 
8 kpc galactocentric radius, thus avoiding the intense star formation in the central kpc. We have recovered 
short-spacing u-v components by using single dish observations from the Nobeyama 45m and IRAM 30m 
telescopes. Using the automated CPROPS algorithm we identified 45 CO cloud complexes in the CO(l — > 0) 
map and 64 GMCs in the CO(2 — > 1) maps. The sizes, line widths, and luminosities of the GMCs are similar 
to values found in other extragalactic studies. We have classified the clouds into on-arm and inter-arm clouds 
based on the stellar mass density traced by the 3.6 /jm map. Clouds located on-arm present in general higher 
star formation rates than clouds located in inter-arm regions. Although the star formation efficiency shows no 
systematic trend with galactocentric radius, some on-arm clouds - which are more luminous and more massive 
compared to inter-arm GMCs - are also forming stars more efficiently than the rest of the identified GMCs. 
We find that these structures appear to be located in two specific regions in the spiral arms. One of them shows 
a strong velocity gradient, suggesting that this region of high star formation efficiency may be the result of gas 
flow convergence. 

Subject headings: galaxies: ISM — stars: formation — ISM: molecules — galaxies: individual (NGC 6946) 



1. INTRODUCTION 

The formation and evolution of giant molecular clouds 
(GMCs), the sites where stars form, are thought to be influ- 
enced by the local properties of the environment and global 
galac t ic structure (IDobbs & Bonnelll 120081: iTasker & Tanl 
2009] lTaskerll20TTt iKrumholz et al.H2009l) . Although the in- 
ternal physical pro perties of molecular c louds have been ex- 
tensively studied dSolomon et al.l IT987t henceforth S87), a 
more detailed understanding of their origin and evolution in 
different types of environments is needed. Maps with resolu- 
tion matched to the scale of individual clouds are the key to 
achieve this goal; however, such maps currently exist only for 
a few galaxies outside the Local Group. Due to the coupled 
evolution between GMCs and their host galaxies, achieving 
a more complete understanding of molecular cloud formation 
and evolution processes represents a key goal for studies of 
galaxy evolution. 

While several mechanisms have been proposed, the pre- 
dominant mode for molecular cloud formation and evolution 
in spiral galaxies is not yet clear. In terms of the arm struc- 
ture, spiral galaxies can be classified as grand-design, most 
of them consisting of two symmetrical spiral arms; multi- 
armed, which show several non-symmetrical arms; or floc- 
culent galaxies with multiple shorter arms. For grand-design 
galaxies, the spiral structure is generated by density waves 
induced by bars or companion galaxies. The influence of 
the density waves in the ISM in these galaxies has been 
commonly mod eled assuming a long l ived, rigidly rotating 
spiral potential dKim & Ostrikerl 120021 IWada & Kodal l2004t 
IDobbs & Bonnelll 120061) . On the other hand, flocculent and 
multi-armed galaxies are thought to develop from local gravi- 
tational instabilities, although flocculent galaxies can present 



weak density waves as well dThornlev & Mundvlll997l) . In 
this case, numerical simulations model the galactic potential 
as time depen dent structures, which can change on a time scale 
of -100 Mvr ( IDobbs & Bonnellll2008l;lWada et aTJl20Tll) . 

The formation of massive molecular clouds in grand-design 
galaxies is attributed to the passing of spiral density waves. 
A particular feature of this type of galaxy is the forma- 
tion of spurs perpendicular to the spiral arms due to the 
shearing of the molecular structures as they leave the arms 
dKim & Ostrikerll2002t IWada & Kodall2004b . In their study 
of the northeastern s piral s egment of IC 342, a grand design 
galaxy. iHirota et ail d201 ll) showed that actively star forming 
GMCs are more massive, have smaller velocity dispersions 
and are more gravitationally bound than quiescent GMCs. 
They interpreted the tendency for the star-forming GMCs to 
lie downstream with respect to the quiescent molecular clouds 
as evidence in favor of cloud a coalescence scenario, where 
massive molecular clouds (~ 10 6 Mq) are formed from co- 
agulation of smaller clouds entering the spiral arms, and the 
star formation is ignited once the excess kinetic energy is dis- 
sipated through inelastic collisions of the cloud s (Hirota et 
al. 20 11). In another grand design galaxy, M51, Koda et al. 
(12009ft suggested that evolution of GMCs is driven by large 
scale galactic dynamics, where the most massive molecular 
complexes (Giant Molecular Associations, GMAs) are cre- 
ated by coagulation of GMCs in spiral arms. As GMAs leave 
the arms, they suffer fragmentation by the strong shear mo- 
tions, which would explain the spur structures observed ex- 
tending from the arms into the inter-arm regions. The remain- 
ing fragmented GMCs are not disassociated into atomic gas, 
and they become part of the molecular cloud population ob- 
served in inter-arm regions which survive until the passage of 
another spiral arm. 
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On the other hand. lDobbs & Bonne ll ( 2008) presented a nu- 
merical simulation of flocculent galaxies. They find that the 
gaseous spiral structure essentially traces the potential min- 
imum, in contrast to the case of grand-design spirals where 
an offset between the spiral shock and the potential mini- 
m um is expected. Accord ing to the simulation s prov ided 
bv lDobbs&Bonnelll d2008b and IClarke & Gittinsl d2006l) the 
most massive structures are thought to be produced in re- 
gions where collision or merging between spiral arms can 
occur, yielding an enhanced star formation in those over- 
dense regions, in contrast to the situation proposed for grand- 
design galaxies where the molecular clouds are formed by 
agglomeration within spiral arms. Conversely, the forma- 
tion of inter-arm molecular gas structures are proposed to be 
the result of local instabilities rather than fragmentation of 
GMA's, and the properties of the clouds are approximately 
similar across the disk, with no offsets between gas and star 
formation expected. High resolution observations of floccu- 
lent galaxies have provided evide nce in support of this picture. 
For instance, iTosaki et al.l (120031) mapped CO( 1 — > 0) over the 
southern arm of the flocculent galaxy NGC 5055. They found 
no obvious offset between Ha and the molecular gas, and that 
on-arm and inter-arm clouds do not have significant differ- 
ences in their properties. 

Resolved GMCs in the Milky Way have been shown to be 
in approximate virial equ ilibrium obey ing scaling relations 
known as Larson's laws (lLarson|[l981l) . These laws estab- 
lish that velocity dispersion (or line width), size, and lumi- 
nosity are correlated in the Milky Way GMCs. The size-line 
width relation may reflect turbulent conditions in the molecu- 
lar ISM. The typical GMC temperature of 15-25 K produces 
thermal CO velocity dispersions of only ~0.1 km s" 1 , provid- 
ing evidence of non-thermal supersonic turbulence within the 
clouds. Nevertheless, whether the turbulence in the molec- 
ular ISM is internally or externally gen erated is still unclear 
dMcKee & Ostrikeril2007t lTaskeril201 ll) . The other two Lar- 
son relations are the luminosity-l ine width and t he lum inosity - 
size relations. In a recent work, iBolatto et al.l d2008l) (hence- 
forth B08) present a study of resolved properties of GMCs 
for a sample of dwarf galaxies along with two disk galaxies 
of the Local Group. They find insignificant differences be- 
tween Milky Way GMCs and GMCs in dwarf galaxies, with 
the latter followin g approximately th e same Larson relations 
as Galactic clouds TBritz et al.l (120071) arrived at a similar con- 
clusion by analyzing observations of CO surveys for galaxies 
in the Local Group. 

Due to sensitivity and resolution limitations, observations 
of resolved properties of GMCs in galaxies other than Milky 
Way are in a n early stage dDonovan Meyer et al1l2012l here- 
after DM12; iKodaet al.l 12009b . Only a few galaxies have 
been mapped in CO emission across the full extent of the 
optical disk with sufficient resolution to probe the ^100 pc 
size s cales on which GMC s form. High resol ution CO sur- 
veys (ISakamotoetal.lll999t iHelfer et al.ll2003l) have focused 
primarily on the central regio ns, whereas comp l ete mapping 
using single-dish telescopes dKuno et al.l 120071: iLerov et all 
120091) has been obtained at resolutions of 11" to 15", corre- 
sponding to ^500 pc for nearby spirals, inadequate to resolve 
GMCs. Only for M33 (BIMA survey) and the Magellanic 
Clouds (NANTEN surveys) has it b een possible to conduct 
a fairly complete census of GMCs (lEngargiolaetal.ll200l 
iFukui et al.ll2008l 

Using a combination of CO(l — > 0) and CO(2 — > 1) and 



TABLE 1 
Properties of NGC 6946. 



Morph. a 
R.A. (J2000) a 
Decl. (J2000 11 
Distance (Mpc) b 
Incl. (°) c 
P.A. (°) c 



SABcd 
20:34:52.3 
60:09:14 
5.5 
33 
243 



Notes. 

a NASA/IPAC Extragalactic Database (NED). 
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pursuing a staged approach of starting with large-scale, low 
resolution maps and following up with high resolution imag- 
ing on smaller areas, CARMA makes it possible to study the 
overall distribution of GMCs, as well as the properties of indi- 
vidual GMCs, in galaxies outside the Local Group. We have 
successfully used this observational strategy in the nearby spi- 
ral galaxy NGC 6946. As indicate d by the single-dish CO 
map from IRAM (ILerov et al.ll2009l) . this galaxy is extremely 
rich in molecular gas, with CO emission extended to a large 
fraction of the optical radius. Table[T]shows the basic param- 
eters of NGC 6946. We have chosen NGC 6946 as our tar- 
get because of its proximity, high CO surface brightness, low 
inclination, and the availability of high-quality datasets at a 
variety of wavelengths, including HI imagin g from THINGS 
and m ultiband Spitzer imaging from SINGS ( K ennicutt et alj 
120031) . The spiral structure in NGC 6946 appears complex, 
with K-band images revealing fou r prominent, asymmetric 
spiral arms dRegan & Vogell 119951) . While the central re- 
gion of the galaxy has been w ell-studied in CO because of 
its high CO surfa ce brightness dCrosthwaite & Turner! 120071 
iKuno et ail 120071 : DM12), the eastern region has not been 
mapped at high resolution. 

In this paper we present observations of the eastern region 
of NGC 6946. Firstly, we made observations of CO(l 0) 
using CARMA covering a region of size 6x6 kpc 2 . The an- 
gular resolution for the CO(l — >• 0) map was ~1 10 pc, which 
is not sufficient to resolve structures with size scales similar 
to GMCs (S87; B08). Then, the second step was to make 
higher resolution CO(2 — > 1) observations (~ 50 pc) towards 
the brightest CO complexes to study the properties of individ- 
ual GMCs. 

We present our study as follows: in Section [2] we describe 
our observations of NGC 6946 using CARMA, and we de- 
scribe the archival data at several wavelengths that we include 
in our analysis. In Section[3]we summarize the technique used 
to identify GMCs and to measure their physical properties. In 
Section|4]we present the scaling relations of the cloud proper- 
ties, and we compare them with previous studies of Galactic 
and extragalactic clouds. In Section[5]we discuss whether the 
properties of the clouds differ between on-arm and inter-arm 
regions. In Section [6] we discuss the implications of our re- 
sults, and in Section [7] we summarize the work presented in 
this paper. 

2. DATA 
2.1. CARMA observations 
2.1.1. 3 mm 

We performed high spatial resolution observations of 
CO(l 0) and 13 CO(l 0) for the eastern part of NGC 
6946 from July to August of 2009. We used the Combined Ar- 
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ray for Research in Millimeter- wave Astronomy (CARMA) in 
E and D array configurations, which have baselines of 8.5-66 
meters, and 11-148 meters, respectively. Our goals are to re- 
solve structures close to the typical sizes of GMCs, and to 
make a detailed study of their properties including sizes, lu- 
minosity, velocity dispersion and star formation activity. In 
order to achieve these goals, we have selected a region in the 
disk with a favorable observing geometry, active star forma- 
tion, and where the interstellar medium transitions from be- 
ing molecular gas dominated to atomic gas dominated. We 
observed CO(l -> 0) and 13 CO(l -> 0) in a 49-point mosaic 
area, which covered 3.6x3.6 arcmin 2 , correspondin g to a 
physi cal scale of 6x6 kpc 2 at a distance of 5.5 Mpc (iTullyl 
1988). The correlator was set to have the CO(l — > 0) line in 
the upper side band (USB), and the 13 CO(l -> 0) in the lower 
side band (LSB). By placing two overlapping 31 MHz bands 
with an offset of 12.9 MHz from the rest frequency line, we 
achieved a total velocity coverage of 150 km s" 1 and a channel 
width of 1.3 km s" 1 . A 500 MHz wide-band was placed out- 
side both spectral windows for calibration. We observed every 
pointing for 30 seconds, yielding a 24.5-minute observation 
time per cycle. At the beginning of the track we observed 
2232+1 17 as the passband calibrator, and at the end of every 
cycle, we observed 2038+513 as the gain calibrator. Calibra- 
tion, imaging and deconvolution were performed using stan- 
dard procedures of the MIR1AD software package. In order 
to have better angular resolution, we applied a Briggs robust- 
ness parameter of -1.5, which yields resolution for CO( 1 — > 0) 
clean map of 5. "2 x 5."0, a cr lms of 0.395 K, and a pixel size of 
1". For the 13 CO(l — > 0) map, because the signal to noise is 
worse than CO(l — > 0), we have applied a Briggs robustness 
parameter of +0.5. The resulting 13 CO(l — > 0) map has a res- 
olution of 5"6 x 5."4, arms of 0.212 K, and a pixel size of 1". 
Figure Q] shows the CO(l — > 0) map o f the region obser ved 
overlaid on a HI image from THINGS (IWalter et al.ll2008l) . 

2.1.2. 1 mm 

Our observations at 3 mm using CARMA likely do not re- 
solve the GMCs. Therefore, the second step was to make 
high resolution observations in D array towards some CO 
complexes to study the properties of individual GMCs using 
CO(2 — > 1). The eleven brightest regions in CO(l — > 0) in- 
tegrated intensity were selected across the area covered by 
3 mm observations. Figure [2] illustrates the regions targeted 
with CO(2 — > 1) observations. The observations were taken 
between July and August 2010. We set the correlator to have 
the CO(2 — > 1) line in the LSB. In this case we use a 125 
MHz band centered on the rest frequency of CO(2 — > 1) line, 
and seven 500 MHz wide bands to observe the continuum at 
1 mm. This correlator configuration yields a velocity cover- 
age of 160 km s" 1 and a channel width of 0.5 km s" 1 . Every 
pointing was observed for 1.5 minutes, which yields a total 
cycle time of 18 minutes, including a 1.5-minute pointing on 
the test calibrator 1927+739. We observed the gain calibrator 
1849+670 after every cycle, and 3C454.3 was used as pass- 
band calibrator. We have applied a +0.5 Briggs robustness 
parameter, which produces a map with angular resolution of 
2 "45x2. "03 corresponding to ~ 50 pc, marginally the scale 
for GMC sizes. The rms noise of the clean maps is 0.527 K. 
The pixel size of the maps is 0.5". The observed regions are 
illustrated in Figure [3] 

We simultaneously obtained data for the same 1 mm re- 
gions using 14 continuum bands covering 3.5 GHz in each 



side band. At 2."43 x 1 ."23 resolution, the rms sensitiv- 
ity is given by a rms = 3 mJybeam" 1 , which corresponds to 
a dust column density Ed us t ~ 7.5 x M Q pc~ 2 , assuming a 
dust temperature of 20 K and a dust opacity at 1.2 mm of 
1 cm 2 g~ l . Further assumi ng a dust-to-gas ratio of 0.01 
(lOssenkopf & Henningll 1994b . this corresponds to a gas sur- 
face density E gas ~ 7.5 x 10 2 M Q pc" 2 . We have no detections 
over the regions observed. 

2.2. Single dish maps 

2.2.1. NRO 45m CO(l -> 0) map 

In order to recover the extended flux, we have merged our 
CO(l ->■ 0) CARMA data with a NGC 6946 single dish map 
from the Nobeyama 45-meter single dish telescope (NRO 
45m). NRO 45m observations of NGC 6946 are part of the 
CARMA and NObeyama Nearby-galaxies (CANON) CO(l- 
0) Survey project, which uses both CARMA and NRO 45m 
telescopes to perform resolved observations of GMCs in disks 
of nearby galaxies (Koda et al., in preparation). DM12 re- 
cently reported observations towards the central part of NGC 
6946 using the CARMA and NRO 45m telescopes, and we re- 
fer the reader to that paper for the observation strategy and de- 
tails of NRO 45m map. Here, we just summarize the most rel- 
evant information. NGC 6946 was observed in the CO(l — > 0) 
transition using the Beam Array Receiver System (BEARS) 
instrument. The FWHM of the 45 meter dish is 15" at the 
CO(l — > 0) rest frequency which is degraded to 19. "7 after re- 
gridding. The velocity resolution of the map is 2.54 km s" 1 . 
The rms noise of the CO(l — > 0) single dish map is 0.13 K 
(0.57 Jy beam" 1 ), which corresponds to ct(Sh 2 ) ~ 5 Mq pc" 2 
assuming a 10 km s" 1 line width window to recover the entire 
flux. 

2.2.2. Heracles CO(2 -+ 1) map 

In the same way as we did for CO(l — > 0) maps, we have 
merged our CO(2 1) CARMA data with a NGC 6946 
single dish map from the HERACLES project dLeroy et alJ 
120091) . HERACLES used the HERA receiver at the IRAM 
30m telescope to map the CO(2 — > 1) transition towards 18 
nearby galaxies over the full optical disk, at 13" angular res- 
olution and 2.6 km s" 1 velocity resolution. At that resolution, 
the cr rms in the data cube is 25 mK, which yields an integrated 
intensity map with a noise level of ct(£h,) ~ 1 Mq pc" 2 as- 
suming that we need to integrate over 10 km s" 1 to recover the 
entire flux. 

2.3. CARMA-Single dish merging procedure 

We have merged our interferometric maps from CARMA 
with the single dish maps from HERACLES and NRO 45m 
using the MIRIAD task immerge. This task linearly merges 
a low resolution image, which represents the short-spacing 
data, with a high resolution image which represents the large 
spacings in the uv plane. The cubes were masked to exclude 
edge velocity channels, and the region where the sensitivity of 
the CARMA map falls below half of the maximum. Immerge 
allows one to solve for the factor to scale the low resolution 
image by to put it on the same flux scale as the high reso- 
lution image. Thus, this factor should be close to 1 if both 
data sets have correctly calibrated flux scales. In our combi- 
nation procedure, we obtained values close to 1.05. We have 
calculated the flux inside the regions of significant emission 
provided by CPROPS (see Section [3] for details) for both the 
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FIG. 1. — CO(l — > 0) integrated intensity contours overlaid on H I integrated intensity map from THINGS. Contours begin at 4tr and are spaced by 4a. The 
noise at the center of the map is 2 K km s _1 . The dashed black line illustrates the region where full gain correction was made. Black dashed ellipses represent, from 
right to left, the region dominated by molecular gas (/ mo i = Eh, /Shi ~ 5; radius = 2.8 kpc), the region where both atomic and molecular gas are comparable 
(/mol ~ 1; radius = 5.6 kpc), and the region where H I starts being the predominant gas phase (/ mo i ~ 0.1; radius = 7.6 kpc). 




RA offset (arcsec) 

FIG. 2. — Left: CO(l — > 0) integrated intensity map. Contours are as in Fi 
observations. Right: 13 CO(l — > 0) integrated intensity map. Contours begin at 1 
Significant emission is observed in regions 5, 6, 7, 9 and 1 1 as is illustrated in the 

image from CARMA observations alone, and the image re- 
sulting from merging CARMA and single dish data. In the 
case of the CO(l — > 0) map, we found that most of the regions 
of significant emission in CARMA+NRO 45m maps present 
< 20% more flux than CARMA maps. Nevertheless there 
are a few regions that can have up to 50% more flux in the 
CARMA+NRO 45m map. On the other hand, we found that, 
for all the regions observed in CO(2 -» 1), CARMA+IRAM 
30m data sets present typically 20% more flux inside the re- 
gion of significant emission. 

2.4. THINGS HI map 

Atomic gas surface density is estimated from 21 cm line 
maps from The HI Nearby Galaxy Survey (THINGS) de- 



RA offset (arcsec) 

gure [T] Black circles indicate the regions where we performed CO(2 — > 1) 
5cr and are spaced by 1.5cr. The noise at the center of the map is 2 K km s . 
left panel. 

scribed in lWalter et ail (120081) . THINGS maps the 21 cm line 
in 34 nearby galaxies using the Very Large Array (VLA). The 
angular resolution of the "robust" weighted map is 4"93 x 
4"51, the channel width is 2.6 km s" 1 , and the noise of the 
map is <~ 17 K. Assuming a 20 km s" 1 integration velocity 
range to fully recover the flux, this noise corresponds to a sur- 
face density of Ehi ~ 7 M©pc" 2 . 

2.5. GALEX FUV data 

NGC 6946 FUV data a re available from the GA LEX Nearby 
Galaxies Survey (NGS; iGil de Paz et ail 120071) . The FUV 
band covers the wavelength range 1350-1750 A, and the 
maps have an angular resolution of 5 "6. We generate a 
sky background subtracted FUV image, by subtracting the 
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sky background image from the intensity map, both provided 
by the NGS project. We examined the image for presence 
of foreground stars over the field we are studying in NGC 
6946. Following the approach presented in lLeroy et al.l (12008b 
(henceforth, L08), we search for foreground stars by using 
NUV/FUV flux ratio. We did not detect the presence of 
foreground stars over the field we are analyzing, so we did 
not apply any masking to the FUV map. We have corrected 
F UV maps for Galactic extinction using the extinction map 
of lSchlegel et al.l (119981) . The FU V extinction estimat ed by 
using A FUV = 8.24 x E(B-V) from lWyder et al.l (120071) . 

2.6. SINGS data 

2.6.1. 3.6 fim map 

We use 3.6 pro, maps fro m the Infrared Array Cam- 
era (IRAC, IPazio et. all 120041) onboard the Spitzer Space 
Telescope to estimat e the stellar surface density (SINGS, 
iKennicutt et al.1 120031) . We aim to identify the spiral arms 
through tracing the old s tellar population. Following a sim- 
ilar approach as taken bv lFovle et al.l (120101) . we have applied 
a spatial filtering to the map to remove variations of 3.6 fim 
surface density which are not related to stellar mass density 
enhancements. We have implemented this procedure by us- 
ing the GIPSY task mfilter, which applies a median filtering 
to the pixels inside a given spatial box. In this study we 
used a 6.75x6.75 arcsec 2 box, which successfully removes 
the bright spots due to recent star formation in the eastern re- 
gion of NGC 6946. 

2.6.2. 24 fxm map 

To estimate the amount of obscured star formation by dust, 
we use the 24 //m data obtained by SIN GS with the MIPS 
instrument of the Spitzer Space Telescope dRieke et al.ll2004l) . 
The angular resolution of MIPS is 5 "7 (FWHM). We perform 
a background estimation as follows. We select a region away 
from the galaxy, and we did a median spatial filtering over this 
background region in order to remove bright spots and other 
features, as we did for the 3.6 fim map. Then, the background 
value is represented by the mean flux over the filtered region. 
The value obtained (~ 0.02 MJy/str) is negligible compared 
to the flux in the eastern region of NGC 6946, thus we did not 
apply a background subtraction to our map. 

2.6.3. Ha map 

In order to estimate the star formation activity of the molec- 
ular structures at the resolution provided by CO(2 — > 1) (~ 
2"), we need a star formation tracer with resolution finer than 
the limiting resolution provided by 24 /im (^6"). With an an- 
gular resolution of ~ 2", Ha maps provide a useful way to es- 
timate star formation for the smallest identified structures. An 
Ha image of NGC 6946 is available in the SINGS fourth data 
release as part of the SINGS ancillary data program. We used 
the stellar continuum-subtracted Ha image, and we convert to 
flux units using the SINGS data release documentation. In or- 
der to remove the [N II] contribution, we have assumed ratios 
of [N II1A6584/Ha= 0.5 and [N II]A6548/[N II]A6584=0.335 
dCalzetti et al.ll2005l) . 

3. CLOUD PROPERTIES 

3.1. Identification 

We identify molecular cloud structures and measure their 
properties using the cloud proper t ies alg orithm (CPROPS) de- 
scribed in iRosolowskv & Lerovl (120061) . CPROPS identifies 



clouds by masking the spectral line cube to isolate regions of 
significant emission that are both spatially and kinematically 
connected, and assigns properties like cloud sizes, luminosi- 
ties and line widths using moments of the emission inside the 
identified region. Those measurements can be corrected for 
the bias due to the finite resolution and sensitivity of the maps. 
This methodolo gy has already been appl ied in several extra- 
galactic studies ?Rosolowskvll2007l; B08: lHughes et al.ll201QL 
iWonget al.1120111) . 

We identify regions of significant emission by selecting pix- 
els with emission greater than a threshold of « t h x er rms across 
two consecutive velocity channels. Then, we extend the re- 
gion of significant emission to the adjacent regions with emis- 
sion greater than « e dge x Ci-ms in two adjacent velocity chan- 
nels. We probed a range of values of «th and « e dge in each 
data set to suppress the presence of noise in our masked re- 
gion. Because the sensitivity in CO(l — > 0) map is better 
than CO(2 — > 1) cubes, we compare the regions identified in 
the latter with the emission found in CO(l — > 0) data cube. 
Regions of marginal significance that do not present a coun- 
terpart in CO(l — > 0) are considered as spurious detections. 
In this study we have used n± = 4 and n^e = 3 for both 
CO(l -> 0) and CO(2 — » 1) maps. 

After regions of significant emission are identified, 
CPROPS will find individual clouds based on conditions that 
depend on the decomposition parameters used. A cloud is 
considered as a separate structure if it is large enough to al- 
low the algorithm to calculate properties, the contrast between 
the peak and the edge is larger than the minimum allowed, 
and the properties of the cloud change significantly when it 
is combined with a nearby identified cloud (see Appendix of 
IRosolowskv & L erov 2006: for a complete description of the 
decomposition parameters). Due to the observed substructure 
in the ISM, the parameters for decomposition must be based 
on the known properties of the objects we want to identify. 
Previous studies hav e used the "physical priors" set of pa- 
rameters presented by IRosolowskv & Lerovl (120061) which are 
based on the observed properties of GMCs in the Milky Way. 
In that case, the decomposition is performed at a spatial reso- 
lution of 15 pc and velocity resolution of 2 km s _1 . Because 
we are limited by the resolution of our maps, in our study we 
use parameters based on the spatial and spectral resolutions of 
the data cube. Then, in finding local peaks in regions of signif- 
icant emission we search in boxes of « 110 pc (5.1") and 2.5 
km s -1 in velocity for the CO(l — > 0) maps and « 50 pc (2") 
and 2.5 km s -1 in the case of CO(2 — > 1). Because the limit- 
ing resolution of the CO( 1^0) maps likely does not allow 
us to resolve individual clouds, particularly the smallest ones, 
henceforth the structures identified in CO( 1 — > 0) maps will 
be referenced as "molecular complexes". On the other hand, 
structures identified in CO(2 — > 1) will be called "clouds". 

3.2. Size, line width and CO flux 

The sizes, line widths and luminosities of the identified 
structures are derived using moments of the brightness distri- 
bution. The rms size o> is determined by the geometric mean 
of the second moments of emission along the major and minor 
axes. The velocity dispersion a v is determined by the second 
moment of the emission along the velocity axis. The flux of 
the cloud is calculated using the zeroth moment of the emis- 
sion in position and velocity. As we noted at the beginning 
of this section, the sensitivity and resolution of the data can 
affect the derived properties of the clouds. For instance, the 
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FIG. 3. — CO(2 — > 1) integrated intensity contours overlaid on the Ha image for all the regions targeted in the eastern part of NGC 6946. The number given in 
the top-right corner corresponds to the region number in the full mosaic as illustrated in the left panel of Figure[2] Contours begin at 5<j and are spaced by 5<r. 
The noise at the center of the map is 2.3 K km s . 
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size of the clouds will be underestimated by the second mo- 
ment of the flux in a map with finite signal to noise. CPROPS 
reduces this bias by extrapolating the size, the velocity width 
and the flux to the case of infinite sensitivity (i.e. brightness 
temperature of the edge of the cloud equal to K). The extrap- 
olation for the size and the veloc ity is linear, while the extrap- 
olation for the flux is quadratic. iRosolowskv & Lerovl (12006b 
recommend correcting for resolution bias as well, which can 
be significant for extragalactic observations. We correct for 
the effect of finite spatial and spectral resolution by subtract- 
ing in quadrature the beam size and the spectral channel pro- 
file from the extrapolated measurements of size and velocity 
width respectively. 

The physical quantities are derived from the moment mea- 
surements following the st andard procedure described in 
IRosolowskv & Lerovl (120061) and B08. Following the defini- 
tion of S87, the cloud size is defined as R = 1.91a r . Although 
this value is ~ 10% lower than the value we would obtain 
from the equivalent area of the cloud (assuming (3 = 1 in the 
density profile p oc r~P), we used this definition in order to 
compare our values to the sizes presented by S87 and B08. 
The CO luminosity Leo is given by 



Leo 



Fco 



K km s 'pc 2 K km s 'arcsec 2 \pc / \ 180 x 3600 
where D is the distance to NGC 6946 in parsecs. 



(1) 



3.3. Mass of the clouds 

Once we have determined the basic properties of the clouds 
(size, line width and luminosity), other cloud properties can 
be estimated by taking combinations of these basic properties. 
For instance, the luminosity -based mass is obtained from Leo 
by using 



Mi u 



Mr; 



= 4.4 



Leo 



K km s-'pc 2 2 x 10 20 cm- 2 (K km s" 1 ) 



i\-i 



(2) 



where Xqo is the assumed CO-to-H2 conversion factor. This 
conversion factor is thought to depend on local ISM prop- 
erties such as met allicity or radiation field dBell et al.ll2006b 
iLerov et ail 1201 lb . In this study we assume Xco = 2 x 
10 20 cm" 2 (K kms" 1 ) -1 , which is broadly used in previous 
extragalactic studies (e. g., B08; L08). In their study of 
GMCs within the central 5 kpc of NGC 6946, DM12 found 
Xqo = 1 .2 x 10 20 cm" 2 (K km s -1 )" 1 within a factor of 2 uncer- 
tainty, which is consistent with the value assumed here. Equa- 
tion |2] includes a factor of 1.36 to account for the presence of 
helium. In this study, we have assumed 1(2 — > 1) = 7(1 — > 0) 
in order to be consistent with the analysis presented in B08. 

Virial masses are normally calculated under the assump- 
tion that molecular clouds are spherically symmetric. How- 
ever, the observed cloud shapes present a more complex 
structure. Gene ralizing the standa r d viri al analysis for 
spherical clouds, iBertoldi & McKeel (11992b calculated the 
virial mass for spheroi d al clo uds. Following Appendix 
A of IBertoldi & McKeel d 19921) the gravitational energy of 
spheroidal clouds is given by 



W 



3 GM 2 

— atao . 

5 R 



(3) 



of the cloud, while a \ and «2 are the density distribution and 
shape-dependent factors respectively. For a p(r) oc r"* 3 power 
law density distribution, ai is given by a\ = (l-/3/3)/(l- 
2/3/5). On the other hand, «2 will change if the cloud presents 
an oblate or prolate form. If the cloud is oblate, an is given 
by arcsine/e, where e = (yl— y 2 ) is the eccentricity of the 
cloud, and y = Z/R is the ratio between the size along and 
perpendicular to the axis of symmetry. If the cloud presents a 
prolate shape instead, an = arcsinh e/e. Ignoring the magnetic 
energy and external pressure terms in the virial equilibrium 
equation, and assuming (3=1, the virial mass is given by 



Mvi 

Mr; 



1040 



km s 



- - 

pc J a2 



(4) 



which corresponds to the expression for the virial mass as- 
suming a spherically symmetric geometry, with a correction 
factor of l/«2- The distribution of an is observed to be narrow. 
We found that the mean of an for CO( 1 — » 0) complexes was 
0.77, with a variance of 0.02 around the mean. On the other 
hand, the mean of «2 for CO(2 — > 1) clouds was 0.81, with a 
variance of 0.02. Thus, the shape factor an can account for a 
correction of the spherically symmetric virial mass by ~ 30% 
for complexes, and ~ 20% for CO(2 — > 1) clouds. 

A complete description of the sizes and ori entations of 
the cl ouds in three dimensions is preferable (e.g. lShetty et alJ 
12010b . but that information is not available for observational 
data. In particular, a correction for inclination can be only 
applied statistically, and not for individual clouds. Thus, we 
assume here that the clouds can be fully described by the sky- 
projected major and minor axes. We make the further assump- 
tion that the axis of symmetry is given by the major axis of 
the clouds, and the clouds are prolate. The correction for the 
shape of the cloud can strongly affect the virial mass estimates 
for structures that present a large axis ratio. Because we as- 
sume elongation in the plane of the sky, a structure elongated 
along the line of sight will appear roughly circular, and our 
correction will be inappropriate. Nevertheless, given the sizes 
found in this regions (<~ 50-150pc), and t he usual thickne ss of 
the molecular gas disk in galaxies (e.g. jYim et al.ll201ll find 
molecular gas thickness ~200 pc for NGC 891), it is likely 
that the axes of symmetry of the clouds are directed along 
the disk, reducing the overcorrection we may be applying to 
our estimates. Supporting this assumption. iKoda et alJ (12006b 
have provided observational evidence that molecular clouds 
are elongated along the Galactic plane. 

We followed the d efinition of the dimension less virial pa- 
rameter a v j r given bv lBertoldi & McKed (fl992l) . 



M lu 



(5) 



where R is the radius perpendicular to the axis of symmetry 



which measures the relative importance of the cloud's kinetic 
energy compared to its gravitational energy. This parameter 
is usually used to investigate whether the clouds are in virial 
equilibrium, a ~ 1 ; for a > 1 the self-gravity is not important, 
and the clouds are confined due to an external pressure if the 
clouds are assumed to be in steady state. 

The uncertainties assessed to the moment measurements 
(and thus to the physical quantities and their derivatives) are 
estimated using bootstrapping methods, and this is the only 
source of uncertainties that we include in our analysis. 

3.4. LTE masses 
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Additionally, we have used the 13 CO(l 0) line map to 
calculate LTE masses over the regions with significant emis- 
sion. If we assume that all energy levels are populated accord- 
ing to local thermodynamic equilibrium (LTE) at the temper- 
ature Tex, the column density of the 13 CO molecule in terms 
of the opacity of the J = 1 — > transition is given by (e. g., 
iRohlfs & Wi£onll200l 



^3 = 3.0X10" 



l-exp(-5.3/T ex ) 



Todv, (6) 



where v is in km s 1 . In this study we have determined the 
excitation temperature by using the assumed optically thick 
line CO(l ^0) 



5.5 



ln(l+5.5/r,i 2 + 0.82)' 



(7) 



and assuming a filling factor of unity. The mean of T ex is 
4.2 K considering the regions of significant emission in the 
CO(l — > 0) map, while this mean increases to 5 K when con- 
siderin g the 1 3 CO(l — > 0) regions of detectable emission. In 
section l4~.2.3l we will investigate the effect of using higher val- 
ues of T ex in the LTE mass estimates. The optical depth of the 
13 CO(l 0) line is given by 



T13 =-ln 



T 13 
1 B 



5.3 Vexp(5.3/r ex )-l 



-0.16 



(8) 



By assuming that 13 CO(l — > 0) line is optically thin, the 
total molecular mass can be computed using an appropriate 
abundance ratio, 



Mi 



Mr. 



3 CO 



D 2 / N( l3 CO)<m 



(9) 



where /i,„ is the mean molecular mass per H2 molecule, D is 
the distance, and [H2/ 13 CO] is the abundance ratio. In this 
study we have assumed the abundance ratio to be 7 x 10 . 
This value is close to the ratio given by the relation found 
by |MilametaII(l2005l) . which yields [ 12 CO/ 13 CO] - 60 for 
^gai ~ 6 kpc, and assuming a H2/" 2 CO abundance ratio of 1.1 
x 10 4 dFrerking et al.|[l982l) . 

3.5. Mass surface density 

Another important cloud property is the mass surface den- 
sity. Surface density has been estimated by dividing the 
luminosity-based mass by the corresponding area covered by 
the CO(2 1) cloud (or CO(l ->■ 0) complex), thus £ H2 = 
M\ um /(irR 2 ). In general the areas for CO(2 — > 1) clouds are 
smaller than CO(l — > 0) complexes. The fluxes and areas in- 
volved in these calculations are the non-extrapolated values. 
We have made this choice, because we intend to compare 
molecular gas surface density with other tracers (HI, 24/im, 
Ha, FUV) over the boundary of the clouds given by the mask 
provided by CPROPS. Because we do not have extrapolated 
or deconvolved measurements of the fluxes for these tracers, 
we have used the raw properties of the clouds. 

3.6. Star Formation surface density 

In this paper, we have estimated the star formation rate 
(SFR) using FUV, 24/im and Ha images described in Section 



[2] The method to estimate the star formation activity is differ- 
ent for CO(l — > 0) complexes and CO(2 — > 1) clouds, mainly 
due to the limiting spatial resolution of the SFR tracer images. 
The details of the calculations are given below. 

3.6.1. SFR in CO(i —¥ 0) complexes 

Following the approach given by L08, we estimate the star 
formation rate surface density in CO( 1 — > 0) complexes using 
GALEX FUV and Spitzer 24 fim maps. By tracing photo- 
spheric emission of O and B stars, FUV emission measures 
unobscured star form ation over time scales of 10-100 Myr 
dCalzetti et al.l 120051) . On the other hand, 24 /j,m traces the 
small dust grains heated by UV photons radiated from young 
stars, thereby tracing obscured star formation over time scales 
of -10 Myr dCalzetti et al.l 12001 iCalzetti et al.1 120071) . The 
use of 2 4 11m as tracer of ob scured star formation was ana- 
lyzed by Calz etti et al.l d2007l) . They showed that the combi- 
nation of Ha (unobscured SFR) and 24/im (obscured) recov- 
ers the total SFR, which they estimated by using Paschen-a 
(Paa) emission as a tracer of ionizing photons unaffected by 
extinction. L08 used a similar prescription to derive the total 
SFR, but using FUV maps instead of Ha to trace unobscured 
SFR. They found that the total star formation surface density 
Esfr can be estimated by: 



^SFR 



M yr _1 kpc 2 



: (8.1 X 10" 2 /fuv + 3.2 X 10" 3 / 24 um)COS! (10) 



where /puv and hA^m are in units of MJy/sr. 

3.6.2. SFR in CO{2 -> 1) clouds 

The coarse resolution of the 24/im map (5 ."7) does not per- 
mit us to use the SFR from Equation dT0t to perform a direct 
comparison between star-forming regions and the structures 
identified in CO(2 — > 1). We therefore use the Ha image to 
estimate the amount of star formation in CO(2 — > 1) clouds. 
However, Ha emission can be strongly affected by dust ex- 
tinction, and a correction should be applied to Ha-based SFR 
estimates to account for embedded star formation. Following 
ICalzetti et alJ (12007ft . the star formation surface density can 
be estimated by 



Ssfr 



M Q yr _1 kpc 



: (5.3 x lO -42 ^) x lO^/^cos/ (11) 



where L Ha is in units of ergs s" 1 , andA Hce is the Ha extinction. 
In the Section 1431 we will show that applying a Au a = 1.0 
mag roughly recovers the total star formation as traced by 
FUV+24/im. 

4. RESULTS 

4. 1 . Uncorrected properties 

Although the CPROPS method corrects for sensitivity and 
resolution bias, those corrections can introduce some addi- 
tional scatter in the cloud property scaling relations. There- 
fore, we will start this section by showing the scaling relations 
without any correction to assess overall correlation strengths. 
We have identified 45 complexes in the CO( 1—^0) channel 
map. The mean uncorrected radius of the CO(l — >• 0) clouds 
is 128.5 ± 3.1 pc, and the mean uncorrected line width is 4.2 
±0.1 km s" 1 . On the other hand, the mean uncorrected lumi- 
nosity of these complexes is (8.6 ±0.1) x 10 5 K km s" 1 pc 2 , 
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TABLE 2 

Spearman rank correlation coefficient r for both raw and 

BIAS CORRECTED RELATIONS. 





CO(l 


— Y 0) complexes 


CO(2 -> 1) clouds 


Relation 


Raw 


Bias Corrected 


Raw Bias Corrected 


a v -R 


0.63 


0.32 


0.60 0.17 




0.77 


0.50 


0.74 0.64 


M m .-L co 


0.87 


0.76 


0.86 0.74 



which corresponds to a mean luminosity-based mass of 3.8 ± 
0.05 xl0 6 M Q . Figure |4] shows the scaling relations for the 
raw properties. 

We selected the 1 1 most massive CO(l — > 0) complexes to 
make follow-up observations with CO(2 — > 1), and they are 
shown in Figure [2] We identified 64 giant molecular clouds 
in the regions observed, and their scaling relations are shown 
in Figure |4] as well. For this set of clouds, the mean radius is 
54.9 ±1.1 pc, the mean line width is 3.7 ±0.1 km s , the 
mean luminosity is (2.4 ± 0.03) x 10 5 K km s" 1 pc 2 , and the 
mean luminosity-based mass is (1.04 ± 0.01) x 10 6 Mq. 

The significance of the correlations between the proper- 
ties are estimated using the Spearman rank correlation co- 
efficient r. We want to assess the change of this coeffi- 
cient when both sensitivity and resolution bias are corrected 
for. Table [2] shows the values found for the scaling relations 
before and after the corrections for resolution and sensitiv- 
ity have been applied. For the CO(l — > 0) complexes and 
CO(2 — > 1) clouds, the <r v -R relation has r ~ 0.6 for the 
uncorrected properties, while the correlation is weaker when 
corrections are applied (r < 0.35). The extra scatter in the 
latter relation can be explained by the resolution bias cor- 
rection. This correction is sensitive to the axis ratio of the 
cloud being deconvolved, and can be extremely large if the 
minor axis is close to the beam size. Thus, elongated clouds 
can deviate more from the original relation than more circu- 
lar clouds of similar average size. For the Leo _ o~ v relation, 
CO(l — > 0) complexes and CO(2 — > 1) clouds show a similar 
correlation r ~ 0.76, and the relations are relatively stable af- 
ter corrections (r ~ 0.65 for CO(2 — > 1) clouds and r ~ 0.5 
for CO(l — > 0) complexes). Among all the scaling relations, 
the M v i r -Lco relation presents the highest r values: r ~ 0.86 
for raw data, and r ~ 0.75 for the corrected values for both 
CO(l -> 0) and CO(2 1). This is expected, as both M vil - 
and Leo involve the product of size and line width. 

4.2. Corrected properties and scaling relations 

Larson's laws reflect the dynamic state of the ISM, which is 
dictated by the interplay between self-gravity, turbulence and 
feedback processes from star formation. Following Larson's 
pioneering work, and using an improved sample of Galactic 
clouds, S87 found that the line widths of the clouds increase 
as a power of their radius, i.e., the size-line width relation is 
given by er v = 0.72 R 5 km s , where er v is the velocity dis- 
persion of the clouds, and R is the radius in parsecs. Com- 
bining this relation with virial equilibrium, they showed that 
all GMCs have approximately the same mass surface density 
£ w 170 M Q pc" 2 . Additionally, they found that the virial 
mass and the CO luminosity follow an approximately lin- 
ear correlation, M v ; r = 39 L C q' Mq, where Leo has units of 
K km s" 1 pc 2 . The third law relates the CO luminosity and 
the line width via Leo = 130 cr, 5 Kkms" 1 pc 2 . 

As observational capabilities have improved in the past 
few years, several studies have established scaling rela- 



tions for extragalac tic clouds as well dBlitz et al.ll2007t B08; 
iHughes et al.ll2.010b . In their study of molecular clouds in Lo- 
cal Group and dwarf galaxies, B08 found scaling relations 
similar to those found in Galactic cloud studies. Neverthe- 
less, they noticed that GMCs in low metallicity dwarf galaxies 
(such as the Small Magellanic Cloud) are larger than clouds 
found in the Milky Way, M31 and M33 for a given velocity 
dispersion or CO luminosity. Th is observational trend was 
confirmed bv lHughes et alj (120101) in their study of the GMCs 
in the Large Magellanic Cloud. 

To find the best fit relations between the properties of the 
clouds, in this pape r we used the bis ector linear regression 
method described in lsobe et al.l (119901) . which consists of the 
line that bisects the ordinary least-square regression (OLS) 
lines in both axes, i.e. OLS(TO) and OLS(XIT) lines. This 
method treats the variables symmetrically, by taking into ac- 
count the uncertainties measured in the independent and de- 
pendent variables, along with the intrinsic scatter of the data. 

4.2.1. Size-line width relation 

Figure|3a) shows the scatter plot of the line width and size 
of both complexes and clouds. Along with the scatter plot, we 
have included the relations found by S87 and B08 for compar- 
ison. The best-fitting relation for CO(2 — > 1) clouds is 



km s 



= (0.14 ±0.05)/?' 



0.88±0.08 



(12) 



where R is in pc. Similarly, for CO(l — > 0) complexes the 
best-fitting relation is 



km s 



■ = (0.27 ±0.19)7? 



i0.65i0.14 



(13) 



We observe that CO(l — > 0) complexes show smaller line 
widths than those found by S87 and B08 for a given radius. 
On the other hand, CO(2 — > 1) clouds are roughly centered 
in the S87 and B08 fitting lines, but present a bigger scatter 
around those relations. Additionally, in Figure|5]we have in- 
cluded the GMCs found by DM12 in their study of the central 
part of NGC 6946 at a spatial resolution similar to our resolu- 
tion forCO(2 — > 1) maps (~ 2"). Although we have identified 
clouds with similar sizes and line widths, they found a set of 
clouds with larger velocity dispersion at sizes ~400 pc. These 
clouds happen to be located within ~400 pc of the center of 
the galaxy. Sim ilar behavior is observed for clouds near the 
Galactic center (lOka et al.ll2001l) . 

4.2.2. Luminosity-line width relation 

Figure |5jb) shows the Lco-ov relation for the clouds in our 
sample. This figure reveals that in NGC 6946, clouds are gen- 
erally located above the relation found by B08, i.e. they are 
overluminous for their velocity dispersion. Nevertheless, this 
trend breaks down for clouds with luminosity smaller than 
~ 3 x 10 5 K km s -1 pc 2 , where the velocity dispersion spreads 
over a wide range. The best-fit relation for CO(2 — > 1) clouds 
yields 



Leo 



K km s pc 



: ((4.99 ±2.99) x io 3 )cc 66±a33 (14) 



where a v is in km s 1 . The best-fit luminosity-line width rela- 
tion for emitting complexes derived from CO(l — > 0) is given 
by 
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FIG. 4. — Properties of the 45 molecular complexes and 64 clouds found in this study without sensitivity or resolution corrections, (a): Size-line width relation. 
Black dots represent the clouds identified in the CO(2 — > 1) maps, while red dots correspond to CO(l — > 0) complexes, (b): Luminosity - line width relation, (c): 
Virial mass - Luminosity relation. Solid blue lines show different Xqq values. 

where Leo is in K km s" 1 pc 2 . On the other hand, the best fit 



Leo 



K km s 1 pc 2 



: ((1.91 ± 1.15) xlOV 



2.41±0.32 



(15) 



Our results for CO( 1 — Y 0) complexes are consistent with 
the values presented by DM12 for the central region of 
NGC 6946. CO(2 1) clouds, on the the other hand, show 
higher velocity dispersions for luminosities below ~ 3 x 
10 5 Kkms" 1 pc 2 . 

4.2.3. Mass estimates of the observed molecular gas 

According to the conceptual model provided by S87, a 
mass-luminosity relation is a natural consequence of the in- 
ternal structure and gravitational equilibrium of the observed 
clouds. The empirical relation found by them (M vlI = 39 L cc \ l ) 
can be explained by combining the size-line width relation 
with the assumption of virial equilibrium, plus assuming sim- 
ilar kinetic temperatures and area filling factors for the clouds. 
This provides a framework for using the optically thick CO 
emission as a tracer of the molecular mass, in which the 
virial masses and CO luminosities of clouds are used to solve 
for the CO-to-H2 conversion factor Xco- Nevertheless, they 
emphasized that this relation depends on the temperature of 
cloud, and may be different in extragalactic sources, particu- 
larly at the centers of galaxies where temperatures and densi- 
ties may deviate from values for molecular clouds in the disk 
of the Milky Way. Surprisingly, recent extragalactic stud- 
ies have found that values of the conversion factor are rela- 
tively close to those derived from Galactic studies. For in- 
stance, for a sample of disk and dwarf galaxies in the Local 
Group, B08 found an Xco consistent within a factor of two 
with the value derived for Galactic clouds. Moreover, using 
high resolution observations of GMCs within the central 5 kpc 
of NGC 6946, DM12 found an X co average value of 1.2 x 
10 20 cm" 2 (K km s" 1 )" 1 , just a factor of 2 below the value de- 
rived from S87. 

The virial mass-luminosity relation for our observations 
is shown in Figure EJc). This figure shows that M v ; r val- 
ues are well predicted by the CO luminosity, supporting 
the use of a relatively constant CO-to-FL factor Xco = 2 x 
10 20 cm" 2 (K km s" 1 )" 1 adopted in equation ©. The best-fit 
relation for CO(2 — > 1) clouds yields 



M v 



^ = (32.2±64.7)L co 



0.87±0.15 



(16) 



relation for CO(l — > 0) emitting complexes is given by 



Mr. 



= (0.12±0.26)L CO 



1.24±0.14 



(17) 



Once again, we have compared the values found in this pa- 
per to the values found by previous works. The structures 
identified in this paper are located in the M v - n - - Leo relation 
found by B08, but are systematically above the relation found 
by S87 for Galactic clouds, and the values found in the center 
of NGC 6946 by DM12. 

In order to further test the feasibility of using CO luminos- 
ity to trace the molecular mass, we have taken advantage of 
our observations of 13 CO(l — > 0). We have detected 6 clouds 
in the 13 CO(l — > 0) intensity map, and their properties are cal- 
culated as we did for CO(l -> 0) and CO(2 ->• 1) clouds. We 
have calculated their LTE masses (Mlte) using the approach 
explained in Section 13.41 The left panel of Figure [6] shows 
the histogram of the CO( 1 — >• 0) complex luminosity -based 
masses, indicating the clouds with 13 CO(l — > 0) counterparts 
detections. We noticed that significant 13 CO(l — > 0) emission 
is detected over several of the most massive CO( 1 — >• 0) com- 
plexes. The right panel of Figure [6] shows the comparison 
between 13 CO(l —>■()) Mlte and luminosity-based mass de- 
rived from CO(l -> 0), both calculated over the 13 CO(l -> 0) 
emission region. We observed that M\ um is usually larger than 
Mlte over the regions we have detections in both CO( 1 — >• 0) 
and 13 CO(l — s- 0) lines. The departure of luminosity-based 
mass from LTE mass is more pronounced for more mas- 
sive clouds. Similar discrepancies between LTE masses and 
other masses tracers have been reported in the literature. For 
instance, in their revision of the Larson's scaling relations, 
iHever etail (120091) found that LTE-derived masses estimated 
from 13 CO(l — > 0) observations of Galactic clouds are signif- 
icantly smaller than the virial masses derived by S87. This 
discrepancy can be explained by the systematic errors yi elded 
by the assumptions in both methods. IHever et all (120091) gen- 
erated models of 12 CO and 13 CO line intensities for several 
physical conditions of the clouds. They found that LTE de- 
rived column densities can both underestimate or overesti- 
mate the true column density, depending on the volume den- 
sity regime, kinetic temperature, among other factors. In ad- 
dition, the assumption of a constant CO abundance inside the 
clouds is another source of uncertainty in calculating the LTE 
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FIG. 5. — Scaling relations for the corrected properties of the complexes and clouds in our sample. After correction, we have 23 CO(l — ¥ 0) complexes and 
25 CO(2 — > 1) clouds in our sample. As in Figuref4] CO(l — ¥ 0) complexes and CO(2 — ¥ 1) clouds are illustrated with black and red circles respectively. The 
black solid line represents the bisector fit for the CO(2 — ¥ 1) clouds, while the red solid line represents the bisector fit for the CO(l — ¥ 0) complexes. The dashed 
line illustrates the relation found by S87, and the dashed dotted line represents the fit found by B08. Red open squares show the values reported by DM12 in 
the central region of NGC 6946. (a): Size-line width relation, (b): Luminosity-line width relation, (c): Virial mass-Luminosity relation. Solid blue lines show 
different Xqq values. 
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FIG. 6. — Left: Histogram of luminosity-based mass for CO(l — > 0) complexes. Black outlines represent the complete sample, while blue shaded regions 
represent the complexes which have 13 CO(l — > 0) detections. Vertical lines illustrate the bin size used in the histogram (10 6 M Q ). 13 CO(l -> 0) is detected 
in several of the most massive CO(l — > 0) complexes. Right: Comparison between LTE masses of 13 CO(l — ¥ 0) clouds and luminosity-based mass of the 
corresponding CO(l — > 0) complexes (black dots). Blue triangles show LTE masses derived using T cx =20 K and red squares illustrate LTE masses derived 
assuming 7^ x =30 K. We observe that M\ um is usually larger (by a factor of ~ 2 or less) than A/lte over the regions we detect in both CO(l -> 0) and 13 CO(l -> 0) 
lines. The departure of luminosity-based masses from LTE masses is more pronounced for more massive clouds. 

- 3 x 10 6 M, 



and luminosity-based masses. The abundance of CO in the 
deeply embedded material of the cloud likely will be different 
from the the material being exposed t o external UV radiation. 
According to Goldsmith et al. (2008) in their study of Taurus 
molecular cloud, the total cloud mass estimated considering 
the CO deficient envelopes can be a factor of two larger than 
the estimated mass assuming a constant CO abundance across 
the cloud structure. Finally, we mention that some difference 
between Mlte and M\ um is expected, since these two differ- 
ent lines, CO(l -> 0) and 13 CO(l -> 0), trace different density 
regimes in the molecular gas. 

Along with the LTE masses from 13 CO(l —> 0) map, Fig- 
ure [6] shows the effect of varying the excitation temperature 
on the derived masses. It shows that it would be necessary 
to increase the excitation temperature to 20 K in order to 
make the LTE mass agree with the luminosity-based mass 
derived from CO(l — > 0) for complexes with masses below 



q. On the other hand, for clouds more massive 
than ~3 x 10 6 M©, a T ex = 30 K is required to account for 
the difference between both mass tracers. Similar T ex has 
been observed in nearb y molecular clouds (Mivawa ki et al.l 
l2009t IWang et ail 120091) . Considering the results presented 
above, we conclude that while the LTE masses are smaller by 
a factor of 2 than luminosity-based masses in our sample of 
clouds, this discrepancy does not mean necessarily that op- 
tically thick CO(l — > 0) observations are not suitable to esti- 
mate true molecular cloud masses, but reflects the systematic 
error related to the assumptions we have made in calculating 
M LTE . 

4.3. Virial Parameter 

Figure [7] shows the virial parameter as a function of M\ um 
and a v for the structures identified in both CO(l — > 0) and 
CO(2 — > 1) maps. The values of a w i r span a range of one order 
of magnitude (0.3 to 2), similar to values that have been ob- 
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FIG. 7. — Left: Virial parameter vs. luminosity-based mass for the iden- 
tified structures in CO(l — ¥ 0) and CO(2 — ¥ 1). Right: Virial parameter vs. 
velocity line width. While a clear correlation between a v j, and M\ am is not 
observed, a v j r shows a correla tion with the veloc ity dispersion. A similar 
correlation has been found by Hirota et al. (2011) in their observations of 
northeastern spiral arm segment of IC 342. Nevertheless, a significant corre- 
lation is expected due to the presence of the velocity dispersion in both axes. 

served in other extragalactic studies (e.g., iHirota et"aDl201 ll 
B08). We observe no clear correlation between a vir and M\ um , 
suggesting that the structur es are bound, rega rdless of how 
massive they are. Although lHirota et al.l (1201 ll) found a sim- 
ilar relation for the virial parameter and the luminosity-based 
mass in their study of the nearby galaxy IC 342, they re- 
ported a good correlation between a v \ r and the velocity disper- 
sion. They suggested that the degree of gravitational binding 
is set mainly by the turbulence within the clouds. Addition- 
ally, they found that clouds associated with H II regions pos- 
sess smaller velocity dispersions and a v ; r than clouds lacking 
H II regions, consistent with a dissipation of turbulence fol- 
lowing passage of the spiral arm. In the right panel of Fig- 
ure |7] we show the dependence of the virial parameter on the 
line width for our observations. Although we observe a weak 
correlation between these two variables, it may be related to 
the fact that these variables are not completely independent 
(avir =Myi r /M\ ma oc <j 2 R/<jTR 2 oc a). 

4.4. Mass surface density 

The mean surface density for CO( 1 — > 0) complexes is 
50 Mq pc" 2 , while for CO(2— >-l) clouds, the mean is 
70. 1 Mq pc" 2 . The difference between the mean surface den- 
sities is not surprising, given that we have selected the bright- 
est regions in the CO(l — > 0) map to observe CO(2 — ¥ 1). 

We note that there is a sub-group of clouds, in both lines, 
that deviate from the main group, which is mostly in the 
regime of E H2 < 100 M Q pc -2 . As we will discuss below, 
those clouds are located in two specific regions in the arms. 
We will discuss the change in cloud properties across the disk 
in Section |5] 

Additionally, we have estimated the atomic gas surface 
density for CO(l — > 0) complexes. Unfortunately, the lim- 
ited resolution of HI maps (~ 5") did not allow us to cal- 
culate the atomic gas component for the CO(2 — >• 1) clouds, 
so we are reporting atomic gas surface densities just for 
CO(l — > 0) complexes. The mean atomic gas surface density 
for CO(l — > 0) complexes is 18 Mq pc" 2 . Figure [8] shows the 
comparison between HI and H2 surface density. We observe a 
maximum value of Shi ~ 40 M© pc" 2 . Moreover, the molecu- 
lar gas fraction, /h2 = Em/Egas, inside the complexes is usu- 
ally high (~ 0.7-0.9), confirming the molecular character for 
most of the regions we are targeting. Nevertheless, there are 
some complexes that show fm below 0.7. Those regions are 
located beyond the radius ~ 4 kpc, where the atomic gas sur- 
face density starts becoming comparable to the molecular gas 
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FIG. 8. — Atomic gas surface density vs. molecular gas surface density for 
CO(l — ¥ 0) complexes. We observe a maximum value of Shi ~ 40 Mq pc~ 2 
(black dashed line) above a molecular surface density of Sh 2 ~ 80 Mq pc -2 . 
The black solid line divides the H I and Hi dominated regions. 

component . Simil ar molecular gas fractions were reported by 
iKoda et ail (12009I) in their study of molecular clouds in M5 1 . 

4.5. Star Formation activity of the clouds 

Having established the properties of the molecular gas 
structures, we now estimate the star formation activity inside 
these regions. GMCs are the nurseries of massive stars in disk 
galaxies, so it is natural to think that the intrinsic physical 
properties of the clouds are tightly coupled to the observed 
rate of star formation. The effect of star formation on the 
prope rties of GMC s has been analyzed through simulations 
(e.g. iTaskerl 1201 ll iMurrav et al.l 120101). t heoretical models 
dOstriker et alj|2010t iKrumholz et al.Jl2010h and observations 
(e.g. IHirota et al .1201 lHHuehes et alJ2010tlWong et al.l20l"lh 
among others. Observations have revealed that star formation 
is a slow process or that a small portion of available gas is 
involved in the formation of stars. Several authors have inter- 
preted this inefficiency in the formation the stars as evidence 
of stellar feedback processes, where the injection of energy 
and momentum into the ISM plays a key role in setting the 
rate at which gas turns into stars. Unfortunately, observations 
of extragalactic molecular clouds usually lack the spatial res- 
olution to relate the star forming regions and their progenitor 
molecular gas, thus making it difficult to interpret the effect 
of star formation in individual clouds. 

Figure [9] shows the comparison between SFR surface den- 
sity for CO( 1 — >• 0) complexes derived from FUV + 24/im and 
Ha + 24/im. As was stated above, the amount of SF is cal- 
culated over the boundary of the identified cloud as yielded 
by CPROPS. We observe that both SFR surface density es- 
timations are well correlated. Along with this comparison, 
we have included the unobscured SFR surface density derived 
from Ha only. The difference between SFR traced by FUV 
and 24/im and the SFR traced by Ha is more pronounced for 
regions of low SFR. In the right panel of Figure [9] we show 
the SFR from Ha corrected for Ah q =1.0 magnitude extinc- 
tion. We notice that this global correction for Ha works rea- 
sonably well in recovering the total SFR (specially for higher 
SFRs), and we have corrected the Equation QT| accordingly. 
This correction is consistent with values found by L08. 
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FIG. 9. — Left: Esfr estimates for CO(l — > 0) complexes as a function of Esf r predicted us i ng a c ombination of FUV + 24/^m, following L08. Black points 
represent the Esfr derived from the combination of Ho and 24/^m as given by Calzetti et al. (2007). Blue points represent the Esfr derived using Ha only. 
Right: Esfr derived from Ha corrected by Ah q = 10 magnitude extinction. 



4.5.1. Esfr vs. Eh2 relation for clouds 

The star formation rate (SFR) is observed to correlate with 
the distribution of gas according to a power law at lar ge scales 
(Esfr ~ E^ as ). This relation was first proposed by ISchmldl 
(119591) . and tested against observations by many authors sub- 
sequently. Th e most influential study of the Schmidt law was 
performed by iKennicuttl (1 1998b . who studied the relation be- 
tween the SFR and gas content averaged over the disk for 
a sample composed of normal spiral and starburst galaxies. 
He found a power-law index N = IA for the entire sample, 
and a steeper index N = 2.47 including just the normal galax- 
ies. As the capabilities of the telescopes have been improved 
over the last years, the study of the Kennicut t-Schmidt law 
has be en performed within individual galaxies. IWong & Blitzl 
(120021) used azimuthal averages of the SFR and surface gas 
density for seven molecular gas rich galaxies, finding indices 
N= 1.2 — 2.1. Using 750 pc resolution d ata for a sample 
of 18 nearby galaxies, IBigiel etaO (120081) found an index 
N=l.O between SFR and molecular surface gas density using 
CO(2 — > 1). They interpreted this linear relation as evidence 
of a constant molecular gas depletion time (~2x 10 9 years) 
within GMCs. Nevertheless, the slope and coefficient ob- 
tained for the Kennicutt-Schmidt law is sensi ti ve to method- 
ology dKennicutt et al.l I2007t iLiu et al.ll20TTt iRahman et all 



in the left panel of Figure [10] We see a clear correlation to- 
wards higher star formation activity for higher molecular gas 
surface density. To find the best fit relation between Esfr and 
Eh2 we used the bisector linear regression method as we did 
for Larson's scaling laws. The best fit relation is given by: 



^ SFR 



M Q yr _I kpc 



:((1.3±0.6)xl<r*)I^ 3±au , 



(18) 



where E h2 is in Mp, pc 1 . This slope is similar to the value 
found by lBigiel et all d2008h (N = 1 .46) for NGC 6946 at 750 
pc resolution for the relation between SFR and the total gas 
component (H I+H2), but is steeper than the relation found for 
H2 (N = 0.92). In Figure[l0]we have included lines of constant 
molecular gas depletion time. Molecular gas complexes are 
consiste nt with a depletion t ime of ~ 1 Gyr, close to the value 
found by lBigiel et all (120081) . 

The right panel of Figure [10] shows the correlation between 
SFR and Eh2 for CO(2 — > 1) clouds. The best-fit relation is 
given by 



^SFR 



M Q yr _I kpc 



= ((1.6±l.l)x 10-=)E 



5svl-96±0.18 



H2 



(19) 



1201 II) . including the fitting method, the spatial resolution and 
the tracers used to estimate the star formation rate and the gas 
density. 

In the present paper, we have chosen to calculate the SFR 
and gas surface density inside the boundaries of the structures 
identified by CPROPS, similar to the approach followed in 
Galactic molecular cloud studies. We aim to establish the 
Esfr-Eh2 relation for scales close to individual molecular 
complex and cloud sizes. We emphasize that this type of anal- 
ysis, spatially biased to the most gas-rich regions, is differ- 
ent from the pixel-by-pixel or azimuthally averaged analyses 
which have been used to d escribe the Kennicutt - Schmidt law 
across galaxy disks (e. g., IWong & Blitzl 120021 : IBigiel et ail 
120081: iBlanc et al.ll2009l) . 

The relationship between SFR surface density and molec- 
ular gas surface density for CO(l — > 0) complexes is shown 



which is steeper than the relation found for CO(l — > 0) com- 
plexes. At these smaller scales, the Esfr-E H 2 relation shows 
higher scatter due to local variations of the SFR and molec- 
ular gas densities. The distribution of points is centered on a 
molecular gas depletion time ~1 Gyr, and shows differences 
of more than 2 orders of magnitude in Esfr between the low- 
est and highest gas surface density extremes. 

4.5.2. Esfr vs. Eh2 using apertures 

By calculating the SFR inside the complexes or clouds, we 
are estimating the star formation activity directly associated 
with these structures. Thus, we may be missing some SFR 
activity nearby to the molecular gas considered in our calcula- 
tions. In order to estimate how the relation between Esfr and 
Eh2 changes when we use larger areas than individual clouds 
or complexes, we have taken an aperture like approach. The 
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Fig. 10. — Left: Star formation rate vs. molecular gas surface density for CO(l — > 0) complexes. The SFR is calculated from the FUV + 24/an map s. The 
black solid line illustrates the best-fit relation found for the complexes, while the black dotted line represents the relation found by Bigiel et al. (2008). The 
vertical dashed line illustrates the sensitivity limit of our map, ~10 Mq pc 2 . Open red squares denote values calculated over 15" (400 pc) apertures centered 
on the molecular complexes. Black dot-dashed lines show constant molecular gas depletion times (SFE~') of 0.1, 1 and 10 Gyr. Right: Star formation rate vs. 
molecular gas surface density for CO(2 — ► 1) clouds. In this case, we have used Ha to estimate the SFR, corrected by 1 mag as is suggested by the right panel of 
Figure|9] Open red squares show the values calculated over apertures centered at the locations of CO(2 — > 1) pointings. The other lines are the same as in the left 
panel, except the sensitivity limit is ~ 15 Mq pc 2 . The fitted slope is steeper than that found for CO(l — > 0) complexes. 



idea is to obtain the enhancement in star formation activity as- 
sociated to the region surrounding the molecular gas clouds or 
complexes, in contra st to th e SFR estimated inside the clouds 
calculated in Section l4.5.1l 

Along with the star formation rate and molecular gas sur- 
face densities calculated inside the structures identified by 
CPROPS, Figure [TO] shows the Ssfr and Eh2 values calcu- 
lated over circular apertures. These apertures have a radius 
of 15" (400 pc), and are centered on the centroids of the 
CO( 1 — > 0) complexes, or at the pointing center of the fields 
observed for CO(2 — > 1) (red squares). We observe that us- 
ing 400 pc apertures, the £sfr-£h2 relation does not change 
significantly with respect to the relation found for individ- 
ual clouds or complexes. In the case of apertures in the 
CO(l — > 0) map, there are a few points with Sh2 ~ 10 Mq 
pc" 2 , but with SFR comparable to other regions with molec- 
ular gas several times larger. Those regions correspond to 
complexes relatively isolated from the main spiral arm struc- 
tures, with moderate star formation activity, and less extended 
molecular gas than SFR tracers. By taking averages over large 
scales, while the molecular gas surface density decreases, the 
Ssfr remains roughly constant. On the other hand, apertures 
calculated over CO(2 — > 1) fields (1 1 in total) roughly follow 
the relation found for individual clouds, although the points 
tend to lie to the left of the best-fitting relation. 

5. ENVIRONMENTAL EFFECT ON MOLECULAR GAS PROPERTIES 

Having determined the properties of the molecular clouds 
in the eastern part of NGC 6946, we turn to the problem of 
investigating if those properties depend on local conditions of 
the gas. More specifically, we search for differences in the 
properties according to the location of the molecular gas. A 
simple way to search for such differences is to use the spi- 
ral structure of disk galaxies. Observations have revealed a 
concentration of young stars in the arm regions of disk galax- 
ies, and given that stars are formed inside GMCs, the prop- 
erties of molecular gas may change once the material enters 



the arm re gion. In their study of the northwestern spiral arm 
of IC 342. iHirota etai] (1201 ll) found that GMCs with active 
star formation are more massive, have smaller velocity dis- 
persions, and are more virialized than GMCs lacking star for- 
mation. Considering that the spatial distribution of GMCs in 
IC 342 indicates that high-SFR molecular clouds are located 
downstream with respect to low-SFR clouds, they suggested 
that GMC properties change following passage through a spi- 
ral arm. In this section, we will take advantage of our spatial 
resolution to compare the properties of the clouds according 
to their location on and off spiral arms. 

5.1. Inter-arm vs. on-arm regions 

Defining spiral arms can be do ne using multiple approa ches 
and techniques. In a recent work. lElmegreen et alJ (1201 lh uti- 
lized arm-inter arm contrast to trace arm structure and the spi- 
ral arm properties of 46 galaxies of different types. They used 
Spitzer 3.6 \im maps due to the small extinction at this fre- 
quency and because the old stellar population dominates the 
emission at this wavelength. Alternatively, spiral arm am- 
plitude can be traced by Fourier transfo rm of the old stellar 
population brightness distributio n ( e.g.. iDurbala et aL l2009t 
IConsidere & AthanassouTal 119881) . lFovleetal.1 ( 120101) used 
Fourier-decomposition of azimuthal scans of deprojected 3.6 
fim maps to define the arm and inter-arm region. In order 
to study how spiral arms affect the properties of GMCs, and 
subsequently the star formation, we have divided the sample 
of clouds into two sub-samples: on-arm and inter-arm clouds. 
We have defined on-arm and inter-arm regions as follows: 
First, we transformed into polar coordinates the deprojected 
3.6 yum image. Table [JJ shows the parameters used to depro- 
ject the images. Then, we decomposed this image into annuli 
with a width of 5 pixels (3 "75). Since the annuli are to be 
Fourier decomposed independently, we overlap the bins by 
70% of a bin width in order to maintain smoothness. For each 
radial bin, the light distribution is expanded in a Fourier series 
by fitting the function: 
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FIG. 1 1. — Identification of arm and inter-arm regions. The image illus- 
trates the 3.6 fim map (color map), and the Fourier-reconstructed mask (black 
contour) defined as enclosing the highest 30% of pixels in each annulus. 



I(r,<t>) = / o (r)+^/ m (r)cos[m(<£-0 m )] (20) 

m=l 

where /„, and </>,„ are the arm amplitude and phase for each 
Fourier component m respectively, and m max is the maximum 
component considered in the expansion. The number of com- 
ponents considered depends on how complex the arm struc- 
ture is. For grand-design ga laxies, spiral arms are well recov- 
ered by using up to m = 4 dFoyle et al.l l2010t> . On the other 
hand, fiocculent galaxies may require expansion out to higher 
components. In order to find the number of components that 
recovers successfully the spiral structure, in our study of NGC 
6946 we have used a range of m max = [5,9]. We observed 
that using ra max = 7 is enough to define the arm structure, and 
adding more components leads to basically the same model. 
Having Fourier-reconstructed the 3.6 /im image, we defined 
the arm region as the area covered by the brightest 30% of 
pixels in each annulus. Figure QT| shows the 3.6 fim image 
(color map) overlaid with the spiral arm mask (black con- 
tour). Once we have determined the spiral arm region, we 
can classify each GMC as on-arm or inter-arm, and search for 
differences in the properties of the two sets. A cloud is de- 
fined as on-arm cloud if the center of the cloud is inside of the 
spiral arm region, and more than the 70% of the cloud's pixels 
are inside the on-arm mask. 

5.1.1. Scaling relations 

Figure [12] shows the histograms of the size, line width, 
luminosity and virial mass of the on-arm and inter-arm 
CO(l — > 0) complexes and CO(2 — > 1) clouds. In the case of 
CO(l — > 0) complexes, we observed that the sizes and line 
widths show no significant differences between on-arm and 
inter-arm structures. We emphasize, however, that due to our 
limiting resolution (~ 1 10 pc at 3 mm) broad line widths ob- 
served in CO(l 0) emitting complexes could be the result 
of high velocity gradient between unresolved clouds. We do 
not observe a significant difference between on-arm and inter- 
arm complexes in the luminosity histogram either, perhaps 
with some on-arm complexes showing slightly larger lumi- 
nosities than inter-arm complexes for a given line width. Fi- 
nally, the CO(l — > 0) virial masses are similar for both sub- 
samples. Similarly, we observe no significant differences be- 
tween on-arm and inter-arm distributions of properties for 
CO(2 — > 1) clouds. Nevertheless, we observe a group of 
clouds located in the high end of all the histograms shown 




FIG. 12. — Histograms of the size, line width, luminosity and virial mass 
for the CO(l — > 0) complexes and CO(2 — > 1) clouds classified as on-arm or 
inter-arm. Black shaded histograms represent structures identified as located 
on-arm, and red shaded histograms show the inter-arm structures. Num- 
bers at the top-right show the number of complexes or clouds associated 
to the on-arm or inter-arm sub-group. We do not observe a clear distinc- 
tion in the properties for on-arm and inter-arm CO(l — > 0) complexes and 
CO(2 — ¥ 1) clouds. Nevertheless, we observe that some CO(2 — > 1) clouds 
show higher luminosities and virial masses than other on-arm clouds (and 
inter-arm clouds). Those clouds spatially match regions of higher star forma- 
tion compared to other locations in the disk. 

for CO(2 — > 1) clouds, as they are more luminous, more mas- 
sive and show broader line widths than the remaining clouds 
located in both on-arm and inter-arm regions. 

5.1.2. Star Formation for on-arm and inter-arm regions 

Having studied the dependence of the properties of the 
clouds on the spiral structure, we will address the ques- 
tion whether the star formation process inside these clouds 
change if they are located at on-arm or inter-arm regions. 
In their study of t hree disk galaxies (including NGC 6946), 
iFovle et ail (120101) investigated the relationship between star 
formation and spiral arms using far-UV and 24/im images as 
star formation tracers. Although they found that the SFR is 
more concentrated in spiral arms, the inter-arm star formation 
activity is still significant (~ 30%). Moreover, they found that 
the fraction of gas in molecular phase is a factor of two larger 
in on arm regions than inter arm regions, although the regions 
of high molecular gas fraction are coincident with regions of 
high total gas surface density in the arms. Thus, this behavior 
of star formation activity and molecular gas fraction provides 
no evidence for a shock-triggered molecular gas formation in 
spiral arms. Nevertheless, the coarse resolution used by Foyle 
et. al (~13") did not allow them to perform a cloud scale study 
of the properties. 

Figure [13] shows the same plot as Figure [10] but using the 
on-arm/inter-arm classification. We observe that on-arm sub 
group is shifted to higher SFR and higher molecular mass sur- 
face density with respect to the inter-arm population. How- 
ever, below Sh2 — 70 M /pc 2 , molecular clouds located in 
inter-arm regions are forming stars at a rate similar to the on- 
arm molecular gas. In order to examine the significance of 
the difference between the SFR surface density of on-arm and 
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FIG. 13. — Ssfr_vs. £h2 relation for CO(l — > 0) complexes (left) and CO(2 — > 1) clouds (right) using the on-arm and inter-arm classification. Lines are the 
same as in Figure [Tol We observe that independent of the spatial resolution, the highest star formation rates and molecular mass surface densities are found in 
on-arm structures. Nevertheless, the typical SFR in inter-arm regions is comparable to that in on-arm regions, especially for low surface density clouds. 

inter-arm clouds, we have used a two-sample Kolmogorov- 
Smirnov (K-S) test. By quantifying the difference between 
the cumulative distributions of two data samples, the K-S test 
allows us to test the null hypothesis that these arrays of data 
values are drawn from the same distribution. A K-S test ap- 
plied to Esfr on-arm and inter-arm distribution of CO( 1—^0) 
complexes yields a D = 0.54 (where D is the maximum dif- 
ference between the cumulative distributions), with a signifi- 
cance of p = 0.0015, while the K-S test for CO(2 1) clouds 
yields D = 0.54 and a significance level of p = 0.0003. Thus, 
because both /^-values are smaller than the default value of the 
level of significance, a = 0.05, we reject the null hypothesis, 
and we conclude that Xsfr for on-arm complexes and clouds 
are significantly different from the inter-arm molecular struc- 
tures. 



5.2. Distribution of the properties in the disk 

In this Section we examine whether the properties of the 
structures identified in this work vary across the disk. Be- 
cause we are not mapping the full extent of the disk of NGC 
6946, we emphasize that the robustness of our analysis may 
be impacted by any large-scale asymmetry in the galaxy. 

A full radial profile study of the ga s and SFR over the opti- 
cal disk of NGC 6946 can be found in lLerov et all ( 120081) . We 
observe that, in general, the size, line width, virial mass and 
luminosity do not show any clear trends as we go from inner 
to outer part of the disk. On the other hand, in Figure [14] we 
present the variation of Eh2, Ssfr and Esfe across the disk. 
In this plot we distinguish between on-arm and inter-arm clas- 
sifications. While there is no overall trend with radius, we ob- 
serve that both the molecular gas and the SFR surface density 
are enhanced for some structures located in on-arm regions. 
This behavior is more pronounced for the SFR, which shows 
clearly a subgroup of clouds that depart from the remaining 
s et of clouds, which shows SFR below ~ 0.06 M yr" 1 pc~ 2 . 

iFoyle et all (120101) found that SFR tracers are more con- 
centrated to the spiral arm region than a uniform distribution 
for the grand design galaxies NGC 628 and NGC 5194, and 
for the galaxy being studied in this paper, NGC 6946. They 
showed that this concentration depends on the pixel fraction 
attributed to arms, and can be 40% using the brightest 30% 



of pixels as the arm region. A similar behavior for SFE is 
observed in the right panel of Figure [14] Several on-arm 
structures are observed to have higher SFE than other struc- 
tures located in on-arm or inter-arm regions. Those regions 
correspond to the highest luminosity clouds identified in the 
histograms for CO(2 — > 1) showed in Figure [12] Spatially, 
those outliers correspond to CO(l — > 0) molecular emitting 
complexes (and the corresponding CO(2 — > 1) resolved inside 
them) located in regions 6, 7 and 1 1 illustrated in Figure[2] In 
order to have a clearer view of the location of these outlier 
structures, in Figure [15] we show the map of the SFE for both 
CO(l — > 0) complexes and CO(2 — > 1) clouds. To create these 
maps, we have assigned the SFE value to the corresponding 
complex or cloud boundary given by CPROPS. The black cir- 
cles highlight the regions where we found the outliers in Fig- 
ures [12] and [14] with black contours showing the CO(l — > 0) 
complexes with Eh2 > 110 M Q pc" 2 , and CO(2 — >• 1) clouds 
with S H 2 > 135 M Q pc" 2 . 

In order to quantify the difference between the structures 
located in regions 6, 7 and 1 1 and the molecular gas in the 
other parts of the regions observed, in Table [3] we show the 
average of Sh2, Ssfr and Ssfe for this sub-group of clouds 
and the remaining set of structures. We observe that the mean 
of Eh2 calculated over the regions located in regions 6, 7 and 
1 1 is a factor of 2 higher than the average value calculated in 
the remaining disk. Also, the average Esfr and Esfe for this 
sub-group are a factor of 4 and 3 higher than the other regions 
respectively. 

6. DISCUSSION 

6.1. Molecular cloud properties 

Despite the long-standing acceptance of the Larson scal- 
ing relations, recent studies have raised questions about the 
form and universa lity of the scaling relations. For instance, 
iHeyer etaf] (120091) . reviewing the properties of the Galactic 
molecular clouds studied by S87, found LTE masses that are 
typically a factor of a few smaller than the virial masses de- 
rived by S87. However, maybe a more remarkable finding 
is that the size-line width relation coefficient (v Q = a v /R 1 ' 2 ) 
is not universal, but depends on the mass surface density of 
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TABLE 3 

Mean properties of the complexes and clouds. 



(Sh2) (Ssfr) 
Mq pc 2 Mq yr~'kpc 



-2 



X 10~ 3 Myr' 



CO( 1 — > 0) complexes in regions 6 ,7 and 1 1 
CO(l — > 0) complexes in the remaining regions 
CO(2 — > 1) clouds in regions 6 ,7 and 1 1 
CO(2 — > 1) clouds in the remaining regions 



121.1 ±13.4 0.17 ±0.02 1.3 ±0.44 

42.9 ± 2.7 0.03 ± 0.002 0.59 ± 0.06 

100.3 ±7.5 0. 19 ± 0.01 2.3 ± 0.43 

57.4 ± 3.1 0.048 ± 0.001 0.79 ± 0.08 



the clouds. iHever etaO (120091) attributed the variation in v D 
to differences in magnetic field strength amo ng the clouds. 
Alternatively, Ballesteros -Paredes et al.l d201 lb proposed that 
the dependence of v on the mass surface density is consis- 
tent with molecular clouds in a process of hierarchical and 
chaotic gravitational collapse. They suggested that, although 
hydrodynamic turbulence is required to induce the first dense 
condensations in the ISM, once the structures become bound, 
the gravity is the main driver of the internal motions. 

Figure [16] shows the relation between v Q and the molec- 
ular gas surface density £h2 for the CO(l — > 0) complexes 
and CO(2 — > 1) clouds. The struct ures we have identi fied ap- 
pear to follow the relation shown in lHever et al.l (12009). where 

1/2 

v oc E H2 is expected for clouds in gravitational equilibrium, 
although these two variables are no independent (S H 2 oc a v ). 

How are the GMCs formed in non-grand design galaxies 
like NGC 6946? In flocculent galaxies (or galaxies with ac- 
tive potentials) the formation of inter-arm molecular gas struc- 
tures is proposed to be the result of local instabilities rather 
than being the fragmented residuals left behind by the kine- 
matic shear of the GMAs in spiral arms of grand design galax- 
ies. In this scheme, the properties of the clouds are approx- 
imately similar across the disk, and no offsets between gas 
and star formation are expected. High resolution observations 
of flocculent galaxies have provided evi d ence in support of 
this scenario. For instance, Tosa ki et"aT] d2003l) mapped the 
CO(l — ► 0) over the southern arm of the flocculent galaxy 
NGC 5055. They found no obvious offset between Ha and 
the molecular gas, and that on-arm and inter-arm clouds do 
not have significant differences in their properties. Similarly, 
in the present study we have identified structures with simi- 
lar properties across the observed portion of the disk. Nev- 
ertheless, we have found that some of the on-arm clouds 
are more massive and have broader line widths than other 
clouds located in both on-arm and inter-arm regions. Why 
are clouds departing from the mean properties observed else- 
where? Accordin g to simulations of ga s in galaxies with 
an act ive potential dClaike & Gittinsll2006l: iDobbs & BonnelH 
l2008tlWada et al.H2011l) the most massive structures are ob- 
served to be produced in regions where collision or merging 
between spiral arms occur, yielding an enhanced star forma- 
tion in those overdense regions. We have found that the more 
massive clouds spatially match sites of higher SFR, which 
may provides evidence in favor of these models. A detailed 
description of the velocity field is needed to examine the po- 
tential prese nce o f colliding flows of gas in this type of galaxy 
(see Section |6T4l below). 

In Figure [12] we found that some of the inter-arm 
CO(l — > 0) complexes can have size and mass comparable 
to the most massive complexes found in on-arm structures. 
However, taking advantage of the finer resolution provided by 
CO(2 — >• 1) observations, we have found that while the most 
massive on-arm clouds remain equally massive as we resolve 
the smaller structures, the inter-arm complexes are decom- 



posed into several less massive components. Nevertheless, 
the sizes, line widths, luminosities and masses are similar for 
most the structures observed in both on-arm and inter-arm re- 
gions. Thus, the properties of the giant molecular clouds in 
NGC 6946 do not change substantially when they enter the 
spiral arms. However, we observe two regions that show sig- 
nificant higher surface density and more massive GMCs than 
in other regions across the disk. 

A natural comparison to our study of resolved GMCs in the 
disk of NGC 6946 may be made with the CO(l -)• 0) obser- 
vations of the central 5 kpc of NGC 6946 presented in DM12. 
As shown in Figure [5^, despite differences in resolution and 
cloud identification algorithms, the sizes and velocity disper- 
sions of the DM12 GMCs and our complexes (as traced by 
CO(l — > 0)) follow a similar trend (with the exception of the 
handful of DM12 clouds within 400 pc of the galaxy cen- 
ter), which is slightly steeper than that measured by S87 for 
Galactic disk GMCs. Figure [5}} shows that the trends de- 
fined by the luminosities of the CO(l — > 0) complexes and 
DM12 GMCs are also consistent. In Figure [5]:, our observa- 
tions indicate that a constant (i.e., not radially varying) value 
of Xqo of 2 x 10 20 is appropriate, which is consistent with the 
value typically assumed for the Milky Way disk clouds. This 
value is slightly higher than (but within the quoted errors of) 
DM12, who derive an average, non-radially varying value of 
1.2 x 10 20 in the central kiloparsecs of NGC 6946. We note 
that adopting an 1(2 —> 1) = 1(1 — > 0) ratio of < 1 would re- 
duce the best-fit value of Xqo for the CO(2 — >• 1) clouds by the 
same factor. Overall, the trend for the CO(l — > 0) detections 
in Figure [5J; is consistent between the two samples, but they 
begin to diverge at the highest cloud masses and luminosi- 
ties, likely because the complexes presented here are blends 
of multiple GMCs. 

6.2. Comparison with Galactic studies 

Recent studies of star formation rate (commonly traced 
by young stellar objects or YSOs) and column densities of 
gas (usually traced by near-IR extinction) in nearby Galactic 
clouds have compared Ga lactic and extragalactic gas and SFR 



surface density relations (|Evans et ai1l2009l : IHeiderm an et al 
l2OT0t iGutermuth et alj|201 ll) . For instance, IHeiderman et al 
(120101) found that the prescription given by the Kennicutt- 
Schmidt law underpredicts the values found in nearby molec- 
ular clouds by factors that range from 21 to 54. They found 
that such differences emerge from the different scales in- 
volved in calculating the SFR and the gas surface densities. 
While the gas tracers used in Galactic studies are usually 
probing the denser component of the gas, the extragalactic ob- 
servations have averaged over scales that include both dense 
and a more diffuse gas component which may not be related 
to the star formation process. 

Figure[l7]compares the values found in th is work for NGC 
6946 GMCs and the values presented by [Heiderman et alj 
d2010l) for Galactic star forming regions. Besides the S gas and 
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FIG. 14. — Radial distributions of £jj2 (left), Sjpr (center) and Esfe (right) of the regions observed over the disk, including both complexes and clouds. 
Filled black dots represent on-arm complexes and clouds, while filled red dots illustrate the on-arm complexes and clouds located in regions 6, 7 and 1 1 (see 
Fi gure s l2l and 1 1 5 1 . Open circles illustrate inter-arm complexes and clouds. Most of the on-arm structures that present higher £h2, ^sfr and Esfe than the other 
complexes or clouds, particularly for the star formation rate, are located in regions 6, 7 and 1 1 . These structures are inside the region enclosed by r = 3.5 kpc and 
r = 4.5 kpc (illustrated by the black dashed lines, see Figure \\5\ . 
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FIG. 15. — Left: Map of the Star Formation Efficiency derived for individual CO(l — > 0) complexes. The color bar is in units of Myr . The black contours 
highlight the complexes with £h2 > 110 Mq pc~ 2 . Right: Map of the SFE for the CO(2 — > 1) clouds. In this case, black contours highlight clouds with 
Eh2 > 135 Mq pc~ 2 . As in the left panel, the color bar is in units of Myr'. Circles illustrate regi ons where we found structures that deviate from the other 
identified structures in Figure [T"4l Dashed lines denote radii of r = 3.5 kpc and r = 4.5 kpc (see Figure fHl . 



Ssfr values for low-mass star formation regions traced by Ay 
maps and counts of YSOs respectiv ely, they in c luded the most 
massive star-forming clumps that IWu et al.l (120101) studied 
with HCN gas maps. We note that low-mass star-forming re- 
gions (red stars) overlap with the GMCs for Eh2 ~ 80 M Q pc 2 . 
Nevertheless, this apparent coincidence has to be considered 
carefully, as we are comparing low-mass (galactic) and mas- 
sive star formation (extra-gala ctic). Additionally, the m asses 
of the large molecular clouds in lHeiderman et alld2010l) study 
are two order of magnitudes smaller than the masses of the 
GMCs reported in this study. Probably a more direct compar- 
ison to our GMCs would be provided by the massive HCN 
clumps, since extragalactic SFR tracers are exclusively sen- 
sitive to massive star formation. In Figure [17] we have in- 
cluded the massive HCN clumps as red diamonds. Surpris- 
ingly, if we extrapolate the £sfr-£h2 re l at i° n we found for 
CO(2 — > 1) GMCs (Equation[l9l) to higher gas surface densi- 
ties, we observe that the massive dense clumps fall roughly 
along the relation. Thus, the £sfr-£h2 re l at i° n f° r massive 
star formation regions is consistent with a quadratic relation. 
Nevertheless, a more complete dynamical range in molecular 



gas and star formation surface density is needed to fill the gap 
between Galactic and extragalactic observations, and assess 
the intrinsic £sfr-£h2 re l at i° n - 

6.3. Star formation and evolution of GMCs in NGC 6946 

Although a complete unbiased survey of the GMCs popula- 
tion and their star formation activity in different environments 
is needed to properly investigate the Xsfr _ relation, our 
analysis allows us to shed light on the difference in star for- 
mati on acti vity betw een on-arm and inter-arm clouds. In Sec- 
tions |3T2] and [372] we observed a clear enhancement of SFR 
in on-arm clouds. This enhancement is more pronounced for 
structures located in regions we suspect are the result of re- 
cent convergence of gas flows (regions 6, 7 and 1 1 in Figure 
[2]). In fact, it is in those regions where we find significant 
13 CO(l — > 0) emission, giving further observational evidence 
for the presence of denser gas in those structures. 

We have found steeper slopes than previous studies on NGC 
6946 for the Ssfr — Sh2 relation for both CO(l — > 0) com- 
plexes and CO(2 — > 1) clouds. As was stated above, the slope 
of the relation Esfr _ can be affected by several fac- 
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FIG. 16. — Variation of the size-line width coefficient v = a v /R l l 2 with the 
mass surface density of the structures derived in this study. The dashed line 
represents the locus of clouds in virial equilibrium. 

tors, including the method used to perform the linear regres- 
sion between these two quantities, the tracer used to estimate 
gas densities and star formation activity, selection effects and 
the resolution of the maps used. For the finest resolution in 
this paper, which is given by the CO(2 — > 1) observations to- 
wards the brightest molecular gas regions, we have found a 
Ssfr- Sh2 relation that is almost quadratic . A quadratic re- 
lation for Esfr-Sh2 has been proposed by iGutermuth et al.l 
d201 ll) for Galactic clouds. They show that their data are con- 
sistent with a star formation law in which the Esfr is pro- 
portional to the E^ as . Deviations from the power law are at- 
tributed to evolutionary stage differences in the local gas: in 
some older regions the gas can be disrupted by stars that have 
been formed there, and other younger regions can contain star 
clusters deeply embedded in dense gas at the onset of star for- 
mation. 

We noticed in Figure Q] that a few regions located in the 
outskirts of the disk present weak or undetectable emission 
in CO(l — > 0) at our sensitivity level, but are bright in 24-fim, 
FUV, Ha and HI. In order to detect the gas in some of those 
regions, we performed deeper low resolution observations 
of CO( 1 — > 0) using CARMA in E array configuration to- 
wards the small region illustrated in Figure [18] The sensitiv- 
ity reached by these CO( 1 — » 0) observations was a ~ 3 M Q 
pc -2 for a channel width of 2.5 km s" 1 and resolution 8. "45 
x 7 "29 (a factor of ~ 4 higher than the sensitivity yielded 
by our extended mosaic). We have applied the same pro- 
cedure explained in Section [3] to identify discrete structures 
and estimate their properties. Nevertheless, in this case we 
are only interested in the £h2, Shi and the Ssfr- Figure [19] 
shows the relation between these quantities for the complexes 
found in this region. The right panel shows that these com- 
plexes exhibit low molecular gas content, but SFR compara- 
ble to denser molecular clouds observed in the complete mo- 
saic (shown by open blue circles). Additionally, in the left 
panel we observe that the atomic content is comparable to the 
molecular gas in these low CO-luminosity complexes. 

A possible way to explain structures with high SFR and low 
molecular gas surface density is to consider an evolu t ionary 
scenario similar to that proposed by IGutermuth et al.l (1201 ll) 
for a sample of large Galactic molecular clouds. Given that 
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FIG. 17. — Comparison between the EsFR - Eh, values found in this work 
for NGC 6946 GMCs from CO(2 -> 1) (black dots) and the values presented 
by Heiderman et al. (2010) for Galactic star forming regions (red symbols). 
Stars show low-mass star formation regions, with Ay maps and counts of 
YSOs used as tracers of gas and SFR r espectively. Dia monds illustrate the 
most massive star-forming clumps from Wu et al. (2010) tr aced by HCN gas 
maps. Red solid lines show the broken power law found bv lHeiderman et all 
120101) . while the black solid line illustrates the relation found in this present 
paper for CO(2 —¥ 1 ) clouds. Galactic low-mass star-forming regions overlap 
with the NGC 6946 GMCs for S H 2 ~ 80 Mq pc 2 , although we emphasize 
that our estimates of SFR trace massive star formation. Thus, HCN massive 
clumps represent a more direct comparison to our GMC values. We observe 
that the massive dense clumps fall along a crude extrapolation of the Esfr- 
Sh2 relation found for NGC 6946 GMCs (Eciuation ll9l to higher gas surface 
densities. 

these molecular complexes in NGC 6946 appear to be iso- 
lated from the spiral arms, they could be GMCs that escaped 
from the gravitational potential minima, and consumed some 
of their molecular gas through forming stars. That atomic gas 
peaks in these regions can be then naturally explained as a 
result of the photo-dissociation of the molecular gas by the 
newborn star population. We observed peak atomic gas sur- 
face densities of ~ 35 M© pc" 2 in these regions, similar to the 
peaks found in the arm regions, and a factor of 3 larger than 
the azimuthally averaged value at that distance from the galac- 
tic center. We postulate that, these structures deviate from the 
main £sfr-£h2 re l at i° n found for the remaining complexes 
or GMCs due to differences in evolutionary stage. Higher 
sensitivity and resolution observations of the gas in these re- 
gions, as well as age estimates are needed to corroborate this 
scenario. 

6.4. Velocity field and comparison with models 

Although a detailed comparison between our observations 
and simulated non-grand design galaxies is beyond the scope 
of this paper, in this section we perform a simple comparison 
between the kinematics of the ISM observed in NGC 6946 
and what is predicted by simulations. In their simulation of a 
galactic disk with a live spiral model. IWada et al] (1201 ll) pre- 
sented a dynamical picture of the behavior of the ISM. They 
showed that both gas and the stellar arms roughly follow the 
local galactic rotational velocity, and they predicted that the 
relative velocities of the gas with respect to the mean rota- 
tional speed should be < 15 kms" 1 . In this scenario, the 
cold gas does not follow the organized flow pattern predicted 
by the conventional density-wave theory and shows rather 
chaotic velocity structures around spiral arm potential min- 
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FIG. 18. — This figure illustrates the region of the eastern part of NGC 
6946 we have observed with higher sensitivity using CARMA in E array. 
The dashed line represent the 7-pointing mosaic we used to cover the region. 
The contour map illustrates the CO( 1 — > 0) intensity map, with contour lev- 
els of 2" X 1.5 K km s" 1 (n > 1). The color map shows the FUV+24/^m 
composite SFR map of the same region. Red ellipses illustrate the molecular 
gas complexes with low CO(l — > 0) emission, but associated SFR. 

ima, especially when the time dependence o f stellar spiral po- 
tential is also ca lculated (IWada et a l. 2011). This effect was 
also observed by lDobbs & Bonnelll (l2008l )in their simulation 
of spiral galaxies with an active potential. 

In order to investigate the kinematics of gas flows in the 
ISM, we have created a velocity residual map in the region 
of NGC 6946 observed in this paper. The velocity residual 
map was created by subtracting the local observed velocity 
field of the molecular gas traced by CO(l — > 0) by a model of 
the global circular motion of the gas. We aim to identify the 
non-circular motions of the molecular gas, and assess whether 
the places with significant velocity gradient are related to re- 
gions of dense gas and high star formation. The model of the 
velocity field was cr eated using the rotation curve derived by 
Ide Blok et all J2008h for the HI map from THINGS. We have 
used the GIPSY task velfi to cr eate the velocity fiel d model 
from the parameters provided bv lde Blok et al.l (120081) such as 
inclination, position angle, and the rotation velocity at a given 
radius using a tilted ring approach. The final residual velocity 
map is shown in Figure [20] We observe that the molecular 
gas shows irregular non-circular motions in the range ~ -10 
km s" 1 to 15 km s -1 in the plane of the disk, in agreement with 
the simulation provided bv lWada et al.l (1201 lb . Moreover, we 
observe a region of strong velocity gradient coincident with 
one of the regions we found to present high SFE (region 1 1 
in Figure ITSb. In order to have a more detailed view of the 
velocity field in this region, in the left panel of Figure I2T1 we 
show a position-velocity diagram of region 1 1 taken paral- 
lel to the local velocity gradient. The molecular gas shows a 
much steeper velocity gradient compared to the global rota- 
tion velocity described by the model. According to models of 
spiral galaxies with an active potential, this region may be the 
result of the convergence of gas into the spiral potential from 
both sides, as we observe in the residual velocity field regions 
moving faster and slower than the global rotation model, i.e., 
a transition of positive to negative velocity values across the 
region. 

The other region with high SFE (region 6 in Figure |2}, on 




FIG. 19. — Left: Atomic gas surface density vs. molecular gas surface den- 
sity for CO(l — > 0) complexes detected with higher sensitivity (red ellipses 
in Figure [T8V We observe that these complexes are located in regions where 
the gas is dominated by the atomic gas phase. Symbols and labels are the 
same as Figure[8] Right: Star formation rate vs. molecular gas surface den- 
sity for the CO(l — > 0) complexes presented in the left panel (black circles). 
As a comparison, blue open circles show th e CO ( 1 — > 0) complexes found 
in the complete mosaic (left panel of Figure [Tol . The comple xes observed 
with h igher sensitivity appear to deviate from the trend found bv. Bigiel et al.l 
12008ft due to their high st ar fo rmation activity but low molecular gas con- 
tent, as is shown in Figure 1181 These cases may represent complexes that 
have been consuming their molecular gas by forming stars, yielding HI from 
photodissociation of Hi by the ambient radiation field. 

the other hand, does not present such a high velocity gradi- 
ent as the region discussed above. However, we observe a 
~15 km s -1 difference in velocity between region 6 and the 
regions located to the north-west and to the south-east. While 
those regions surrounding the molecular gas complex in re- 
gion 6 are moving more or less at -5 km s" 1 with respect to 
the rotating frame, region 6 is moving in the other direction 
with a velocity of ^10 km s . This difference in velocity can 
be observed clearly in the position-velocity diagram shown in 
the right panel of Figure [2]] This may imply that the gas ac- 
tively forming stars in region 6 could be the result of a recent 
convergence of flows. Although we do not observe signatures 
of a merging or collision process of spiral arms in the gas 
distribution near region 6, it is possible that the process had 
already formed a new spiral arm, obscuring the original arm 
locations. More detailed analysis of the kinematics of the gas 
is required to disentangle the actual process involved in the 
enhancement of star formation efficiency in region 6. 

6.5. Caveats of this work 

The study presented here can be improved in multiple ways. 
Our analysis has been limited to the north-eastern area of the 
molecular gas disk of NGC 6946 for CO( 1—5-0) observations, 
and to some of the brightest regions for CO(2 —5-1) follow-up 
observations. A full census of the GMCs population across 
the disk is needed in order to perform a detailed study of 
molecular cloud properties and their relation with the sur- 
rounding environment. Nevertheless, high sensitivity maps 
with resolutions ~ 60 pc remain observationally challenging 
to obtain, although recent studies have achieved this goal with 
CARMA+NRO 45m (e.g. DM12). New facilities like ALMA 
will provide maps of unprecedented quality, allowing us to 
perform an unbiased analysis of the molecular gas in nearby 
galaxies at scales close to GMC sizes. 

In our SFR calculations for CO(2 — > 1) clouds, we have ap- 
plied a crude global extinction correction Ah« = 1 mag to the 
SFR(Ha). Although we showed that such correction is a good 
approximation to the SFR values found using a combination 
of FUV+24/im for CO(l — » 0) complexes, a more sophisti- 
cated approach to calculate the obscured SFR at resolution 
close to the GMCs sizes is needed. Thus, tracers relatively 
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FIG. 20. — Residual velocity field of the region observed in NGC 6946. 
The color bar is in units of km s _1 . Black solid lines illustr ate t he slices 
used to generate the position-velocity diagrams shown in Figure l271 We have 
subtracted a model of the rotation velocity from the observed CO(l — > 0) ve- 
locity field. The model of the rotation velocity was created using the rotation 
curve derived by de Blok et al. (2008). We observe that there is a significant 
velocity gradient in one of the northern regions where we detect high SFE 
(region 1 1 in Figure[l5). In the southern region (region 6 and 7 in Figure [T31 . 
on the other hand, we do not observe a significant velocity gradient inside the 
individual CO complexes, but there is a ~ 15 km s _I discontinuity between 
the gas in regions 7 and 6. 

unaffected by dust absorption represent a valuable alternative 
to estimate the total SFR at smaller scales. In this direction, 
observations of Paschen a (Pa) line could provide unbiased 
measurements o f the number of ioniz ing photons in these star 
forming regions (Calzetti et al. 20Q3). 

Equations \TU\ and QT| assume that stars have been formed 
continuously. Nevertheless, this assumption may not be 
longer valid for scales close to GMC sizes, as one approaches 
the case of a single stellar population with a discrete age. 
Some studies have investigated the uncertainty in using a 
linear conversion from UV or Ha luminosity to estimate 
the SFR, making use of models assuming an instantaneou s 
burst of star format ion dGenzel et al.l2010tlLerov et al.l2012h . 
iLerov et all (120121) showed that while Ha emission occurs 
mostly in the first 10 Myr, FUV presents low, but significant, 
emission up to 65 Myr after the burst. They report a factor 
of 2 uncertainty inferring the SFR from Ha, and a factor 3- 
4 uncertainty inferring SFR from FUV. Resolved star cluster 
observations in the regions observed in this study could shed 
light on the ages involved, and provide an independent way to 
estimate the star formation rate which can be compared to the 
SFR traced by FUV or Ha. 

7. SUMMARY 

We have performed one of the most detailed observations 
of GMCs outside the Milky Way. Using CARMA, we have 
observed CO(l — > 0) and 13 CO(l — > 0) over the northeastern 
disk of the non-grand design galaxy NGC 6946, covering a re- 
gion of 6x6 kpc 2 at a resolution of ~1 10 pc. Although at this 
resolution it is possible to resolve some of the largest molecu- 
lar clouds, we have complemented this map with higher (^50 
pc) resolution maps of CO(2 — >• 1) towards the brightest re- 
gions observed in CO(l — > 0). The extended flux has been re- 
covered by merging our interferometric maps from CARMA 
with single dish maps from NRO 45m and IRAM 30m. The 
results are summarized as follows: 

1 . We have estimated the properties of the molecular gas 





FIG. 21 . — Position- velocity diagrams of the regions showing velocity fea- 
tures in the residual velocity field. The exte nt and orientation of the slices are 
shown by the black solid lines in Figure|20] Color map shows the CO( 1 — > 0) 
data, and the black contours ill ustrate the model of the rotation velocity cre- 
ated using the rotation curve by de Blok et al. (2008). Left: Position- velocity 
diagram of region 1 1 . We observe that the velocity gradient is sharper in the 
CO(l — > 0) data than the global rotation velocity model. Right: Position- 
velocity diagrams of the regions 6 (peak on right) and 7. We do not observe 
a significant velocity gradient inside the individual CO complexes, but the 
molecular gas in region 6 is moving about ~ 15 km s -1 away from the con- 
tiguous region 7. 

structures (size, line width and flux) using the cloud- 
finding algorithm CPROPS (Section g}. We have re- 
solved the largest molecular complexes into several 
structures. At the resolution offered by CO(l — > 0) ob- 
servations, we have identified 45 CO emitting com- 
plexes with typical sizes of ^150 pc. Five of these com- 
plexes show significant 13 CO(l — > 0) emission, which 
appear to be regions of denser gas. The higher reso- 
lution observations of CO(2 — > 1) towards the brightest 
regions detected in CO( 1 — >• 0) allowed us to find 64 
structures which we have identified as GMCs. 



In Section 14.21 we observed that the properties of the 
clouds and complexes follow relations similar to those 
found previously for extragalactic clouds (Figure|5]l. In 
particular, the trends defined by the CO(l — > 0) com- 
plexes presented in this paper are consistent with those 
found in the center of NGC 6946 by DM12, despite the 
differences in resolution and cloud identification algo- 
rithms; the exception is the set of GMCs located within 
400 pc of the galactic center with large velocity dis- 
persions (> 10 km s -1 ). Additionally, our M v j r -Lco re- 
lation is located slightly above the relation found by 
DM12 for the central part of NGC 6946 and is consis- 
tent with our choice of the CO to H2 conversion factor 
X co =2 x 10 20 cm- 2 (Kkm s" 1 )" 1 . 

We have estimated LTE masses for the regions with sig- 
nificant 13 CO(l — > 0) emission (Figure |6). Although 
the LTE masses are smaller by a factor of 2 than the 
luminosity-based masses in our sample of clouds, this 
discrepancy does not necessarily mean that CO(l — > 0) 
intensities are not suitable to estimate true molecular 
cloud masses, but may reflect systematic errors arising 
from our lack of knowledge of the excitation temper- 
ature and CO abundance, both needed to infer masses 
from 13 CO(l — > 0) emission. 

The star formation rates o f the CO(l — > 0) complexes 
were estimated in Section 1431 using multiband Spitzer 
imaging from SINGS and FUV GALEX maps of NGC 
6946. In Figure [10] we observe a clear correlation 
between SFR surface density (Ssfr) and molecular 
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gas surface density (Em). The bisector linear regres- 
sion for the relation Esfr = A yields a slope of 
N =1.43 ±0.11 fortheCO(l -> 0) emitting complexes, 
which is st eeper than the valu e found by previous work 
(N = 0.92, iBigiel et al.ll2008h . In order to detect faint 
CO complexes, but with active star formation, we have 
observed with higher sensitivity the CO( 1^0) line 
over a smaller region located beyond the spiral arms, 
where atomic gas is comparable (or dominant) to the 
molecular gas (Figure [TSll. We interpret these CO-faint 
star forming regions as isolated structures that could 
have departed from the main spiral arms, and have been 
turning their molecular gas into stars since then. Thus, 
it is not surprising to find atomic gas peaks in these re- 
gions, which could have been produced by photodisso- 
ciated molecular gas. Higher sensitivity observations 
are necessary to detect faint molecular gas associated 
with such late stage star formation. 

5. For GMCs identified in the CO(2 — > 1) maps, we es- 
timate their SFR using Ha maps. We have found a 
steeper slope than the values found for CO(l — > 0) com- 
plexes, with N = 1.96 ±0.18. We observe that Galactic 
massive dense clumps follow a rough extrapolation of 
the S S fr-Sh2 relation for these CO(2 1) GMCs to 
higher gas surface densities (Figure |T7T>. ALMA will 
observe the molecular gas distribution at even higher 
resolution and sensitivity, allowing a direct comparison 
between the Galactic and extragalactic £sfr-£h2 rela- 
tions at similar scales. 

6. In Section [57X1 on-arm and inter-arm regions have been 
defined based on stellar mass density enhancements 
traced by the 3.6 fim image of NGC 6946 (Figure QTJi. 
The complexes and GMCs were classified as on-arm 
and inter-arm structures accordingly. Although we ob- 
serve that cloud properties have similar values indepen- 
dent of location in the disk, a small number of clouds 
are more massive, present broader line widths, and 
higher star formation rate than the remaining sample 
of clouds (Figure [12} , We found that these structures 
appear to be located at two specific regions in the spiral 
arms (region 6-7 and region 1 1 in Figure[2]i. 

7. The star formation efficiency (SFE, defined as the ratio 
between Ssfr and J^m) is roughly flat over the region 
observed , wi th values between 10~ 4 and 10~ 3 Myr -1 
(Section 15.21 Figure [141 . However, the SFE is en- 
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hanced in the regions where we observed more massive 
and more luminous GMCs (see FigurefTBIand TableO. 
The velocity field of the molecular gas reveals irreg- 
ular non-circular motion with values below 15 km s -1 
with respect to the rotating frame (Figure l20ll. In ad- 
dition, we found some evidence for convergence of gas 
flows in one of the regions of high SFE (Figure[2TJi. The 
other region presenting SFE peaks does not show such 
a strong velocity gradient. However, this region shows 
a 15 km s -1 difference with respect to the surround- 
ing molecular gas, which could be caused by a previ- 
ous convergence of gas flow. Although a more detailed 
comparison with simulations is needed in order to have 
a more conclusive picture of the role of gas dynamics in 
the formation of GMCs and stars in these type of galax- 
ies, these new observations suggest that, in agreement 
with numerical simulations of non-grand design spiral 
galaxies, the convergence of gas flows into the spiral 
arm, and collisions or merging of arms may increase 
the local gas density to higher values than the mean 
observed in the disk, creating more massive molecular 
clouds, and enhancing the rate of star formation. 
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